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FOREWORD 


This document was prepared by Rocketdyne, a Division 
of Rockwell International Corporation, in accordance 
with DRLT-1204, Line Item No. 4 of the Data Requirements 
List of Contract NAS9-14801 with the National Aeronautics 
and Space Administration. The contract period of per- 
formance was December 1975 to July 1977. The NASA/ JSC 
Technical Monitor was Mr, Ronald C. Pasadyn. The 
Rocketdyne Program Manager was Mr. R, D. Paster for 
the first ten months; he was replaced by Mr. H. S. Stratton 
for the remainder of the program. Mr. J. A. Nestlerode 
served as the Project Engineer and the Principal Engineer 
was Dr, R. C. Kesselring. Mr. R. L. Nelson was responsible 
for the model formulation and programming. Mr, K. M. Sprouse 
was responsible for the analytical descriptions of the in- 
jector modelling technique and the chamber pressure input. 

This report has been assigned Rocketdyne Report No. R77-159. 




ABSTRACT 


An analytical injector model was developed specifically to analyze combustion 
instability coupling between the injector hydraulics and the combustion pro- 
cess. This digital computer dynamic injector model will, for any imposed 
chamber or inlet pressure profile with a frequency ranging from 100 to 3000 Hz 
(minimum) accurately predict/calculate the i'nstantaneous injector flowrates. 

The engine hydraulic stabiity computer model was developed using the "lumped 
parameter" technique. The model was formulated such that the equations govern- 
ing fluid flow inside the injector were linearized and subsequently arranged in 
matrix form. The matrix is then solved in a frequency response format giving 
gains and phases for pressures and flowrates at various locations within the 
injector. 

The injector system is described in terms of which flow segments enter and 
leave each pressure node. For each flow segment, a resistance, line lengths, 
and areas are required as inputs (the line lengths and areas are used in 
determining inertance) . For each pressure node, volume and acoustic velocity 
are required as inputs (volume and acoustic velocity determine capacitance) . 

The geometric criteria for determining inertances of flow segments and capa- 
citances of pressure nodes was set. Also, a technique was developed for 
analytically determining time averaged steady-state pressure drops and flow- 
rates for every flow segment in an injector when such data is not known. 

These pressure drops and flowrates are then used in determining the linearized 
flow resistance for each line segment of flow. 

Model output includes both absolute and vector summations of the oscillatory 
injector flow (for a unit pressure perturbation) so as to facilitate interpre- 
tation of the model results. Graphical displays (plots) are also included. 

The model was correlated with experimental data from three injectors - the 
Rocketdyne Lance XRL Booster injector, the Aerojet Space Shuttle OME 
Technology injector, and the Rocketdyne Space Shuttle OME Technology injector. 
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Each of these injectors experimentally exhibited hydraulic coupling in the 
100- to 3000-Hz range. In addition, experimentally proven hardware fixes 
which successfully eliminated the instability, exist for each of the three 
"correlation" injectors. 

The engine hydraulic stability computer model was run for each of the "cor- 
relatipn" injectors with both prei-fix and post-fix injector model input. 
Examination and analysis of the model output revealed that the computer model 
successfully predicted that the fixes applied to each correlation injector 
would increase combustor stability with respect to the instability mode 
actually observed. 

The establishment of injector design criteria was attempted by conducting 
a sensitivity analysis with the model through a systematic study of various 
injector design variables using the model to constantly gauge the effects 
relative to injector gain (stability). 

It is concluded that the engine hydraulic stability computer model is quite 
capable of predicting effects of proposed injector fixes and/or evaluating 
design options relative to one another. 

The program was completed with the successful operation of the engineering 
model on the NASA/ JSC Univac 1110, EXEC-8 computer system and by extensive 
documentation of the model in the form of a computer manual and final 
report. 
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SECTION I ' 


INTRODUCTION AND SUMMARY 

Combustion instability normally falls into one of three types: acoustic, feed- 

system coupled, or hybrid. 

These types of instability have been classified according to the wavelength of 
the frequency compared to a characteristic dimension of the chamber. When the 
frequency is so low that pressure fluctuations are felt (in bulk) simultane- 
ously throughout the chamber, the instability is directly related to the feed 
system effects and is commonly called "feed-system coupled" or "chug." When 
frequencies become great enough that wave motion within the chamber is involved, 
the instability is classified "acoustic" if feed system effects are minimal, 
and "hybrid or combined" if feed system effects are predominant. The most clas- 
sic instabilities are chug and acoustic, but recently hybrid instabilities have 
been observed in a number of engines. 

Generalized models dealing with acoustic and feed-system coupled instability 
are currently in use. Such a generalized model for hybrid instability, however, 
does not currently exist. Nonetheless, limited models developed for specific 
cases of hybrid instability have achieved success in the past. This program 
is concerned with the development of a generalized model whose emphasis is on 
the hybrid type of combustion instability. 

Analytical models of acoustic instability describe the variation of burning 
rate with position in the combustion chamber. This is shown schematically in 
Fig. 1 where the portion of the engine system being considered analytically is 
enclosed by the heavy black line. The analytical model of acoustic instability 
employed most often at Rocketdyne is a modification of the Priem model. 

As indicated in Fig. 2, the essence of the feed system model is a detailed one- 
dimensional treatment of the feed system and a simplified treatment of the 
combustion process and chamber dynamics. Rocketdyne recently completed such a 
generalized model under contract NAS9-14515 with NASA-JSC (Ref. 1). The feed- 
system model predicts relative stability and reveals the essential design 
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Figure 1. Acoustic Instability - Variations in Burning Rate 

characteristics influencing that stability. An example of this is shown in 
Fig. 3. Fig. 3 is a Bode plot and is widely used for studying system stability. 
The y-coordinate, w/2 AP, is termed the gain, while the x-coordinate, f^, is 
termed the break frequency. Lower gain systems are most stable (other factors 
equal). The effect of AP in lowering the gain is shown. Also shown is the 
effect of increasing orifice inertia or increased orifice L/D. Increasing the 
orifice L/D decreases the break frequency and enlarges the region of low gain. 

The resonant modes are determined by the overall feed system design details. 

The complexity of hybrid instability is indicated in Fig. 4. Two- or three- 
dimensional aspects of both the feed system and the chamber acoustics must be 
taken into account. This is because constituent portions of the injector such 
as ring grooves or annular manifolds have characteristic frequencies which are 
close to the acoustic frequencies of the chamber itself. While this total model 
has not yet been solved analytically, engines exhibiting these phenomena have 
been analyzed and their problems resolved by limiting the analysis to the feed 
system, combustion process, chamber dynamics, or a simplified combination of the 
above . 
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Figure 2. Feed System Coupled Oscillations-Variations in Injected Flowrate 
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Figure 3. Evaluation of Feed System Modifications 








Figure 4, Hybrid Instability 

The digital computer injector model developed during the subsequently described 
effort analyzes an injector in terms of its coupling potential with the chamber 
dynamics after the injector’s multidimensional description has been input to the 
model. The model is an open- loop type model requiring that a combustion chamber 
pressure profile representing a specific acoustic mode of instability be imposed 
on the injector. The model predicts the injector response in terms of flowrate 
and pressure oscillations throughout the injector, and is capable of analyzing 
frequencies up to at least 3000 Hz. 

This document is the final report of a 16-month analytical research program to 
develop generalized injector design criteria for the prevention of coupling be- 
tween the injector hydraulics and the combustion process. This program was 
conducted in the following manner: 
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• Construction of a digital computer dynamic injector model which, 
for any imposed chamber or inlet pressure profile with a frequency 
ranging from 100 to 3000 Hz, accurately predicts/calculates the 
instantaneous injector flowrates, 

• Correlation of the model with existing test cases which exhibited 

hydraulic coupling in the 100- to 3000-Hz range and verification that 
the model can predict characteristic differences in the frequency 
response of injectors: which coupled and those which did not couple. 

• A sensitivity analysis which consisted of a systematic variation of 
injector design variables using the model whose end objective was the 
establishment of generalized design criteria. 

The analytical and test evaluation range of interest in this program is summar- 
ized below: 

• Coupling problem in 100 to 3000 Hz frequency range 

• Propellants - acid/amine. LOX/amine, acid/hydrocarbon, 

LOX/hydrocarbon 

• Thrust - 25 to 50,000 lbf 

• P - 100 to 1000 psia 

• HR - maximum performance +20% 

• Fuel temperature - 40°F to 0° subcooling at injected conditions 

• Acid temperature - 40°F to 0° subcooling at injected conditions 

• LOX temperature - -298°F to 0° subcooling at injected conditions 

• Orifice size - 0,020 to 0.040 in. diameter 

• Orifice AP - 25 to 50 psid (or 15% of P ) 

The specific end products of this effort are as follows: 

1. A practical digital computer program (model) running on the 

NASA/ JSC computer which describes the Injector manifolding and 
orifice feed geometry and permits analysis of the injectors, 
hydraulic coupling characteristics, The model will be structured 
in a manner compatible with a total system stability model (NAS9-14315) . 
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The necessary design criteria which will allow the injector 
designer to configure an injector such that hydraulic coupling 
will not occur. 

3, Complete documentation consisting of the following: 

a, A final report describing the entire effort. 

b, A computer manual with appropriate sections for the 
user, the engineer, and the programmer. 

The program was accomplished in four tasks : 

Task I. Review Existing Models 
Task II, Model Formulation 
Task III, Model Correlation 
Task IV. Sensitivity Analysis 

In Task I - Review Existing Models, data was compiled on the modeling tech- 
niques employed in existing injector models. Three basic injector modeling 
techniques (lumped parameter, continuous parameter, and multidimensional 
wave solution, i.e,. Green's functions) were evaluated for possible utiliza- 
tion in the generalized injector model to be developed in this program. To 
assist in the evaluation, a total of three injectors (all of which have 
exhibited hydraulic coupled instability) were defined in detail. The lumped 
parameter injector modeling technique was demonstrated to be capable of 
satisfactorily describing any of these "typical" injectors and was selected 
for utilization in Task II. Application of this technique to "typical" 
injector designs is simple, flexible, and state-of-the-art. Effort conducted 
in Task I resulted in the elimination of possible disadvantages previously 
believed associated with the lumped parameter technique (i.e., (1) computer 
core storage, (2) computation time, (3) limitation on injector description 
capability due to core storage limitation, and (4) accuracy due to descrip- 
tion limitation. 

In Task II - Model Formulation, the model was formulated such that the 
equations governing fluid flow inside the injector were linearized and 
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subsequently arranged in matrix form. The matrix is then solved in a 
frequency response format giving gains and phases for pressures and flow- 
rates at various locations within the injector. 

The injector system is described (using the lumped parameter technique) in 
terms of which flow segments enter and leave each pressure node. For each 
flow segment, a resistance, line lengths, and areas are required as inputs 
(the line lengths and areas are used in determining inertance) . For each 
pressure node, volume and acoustic velocity are required as inputs (volume 
and acoustic velocity determine capacitance) . The geometric criteria for 
determining inertances of flow segments and capacitances of pressure nodes 
was set. Also, a technique- was developed for analytically determining time 
averaged steady-state pressure drops and flowrates for every flow segment 
in an injector when such data is not known. These pressure drops and flow- 
rates are then used in determining the linearized flow resistance for each 
line segment of flow. The model was constructed to have an input option of 
any of three chamber instability mode shapes (first or second tangential or first 
radial) . However, any general pressure profile input can be used by specifying 
the amplitude and sign of each of the individual input locations. 

To provide a simple way to evaluate any specific injection configuration, 
summations were made of the oscillatory injector flow (for a unit pressure 
perturbation). These summations include: (1) a summation of all the 

absolute values of injector flowrate, (2) a vector summation of all injector 
flowrates (attempt to include phase angle), and (3) and (4) the previous 
two cases except that each injector flow is multiplied by the relative 
amplitude of the chamber pressure that it feeds before making the summations. 

In addition to printout of the summation of injector flows (both absolute 
and vector), the ease of interpreting the output of the engine hydraulic 
stability computer model was also enhanced by a graphical display of the 
amplitudes of ring groove flows, ring groove pressures, injector flows, and 
chamber pressure inputs around each ring groove. 
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In Task III - Model Correlation, three "correlation" injectors were selected 
for utilization. These were the Rocketdyne Lance XRL Booster Injector, the 
Aerojet Space Shuttle OME technology injector, and the Rocketdyne Space Shuttle 
OME technology injector. Each of these injectors experimentally exhibited 
hydraulic coupling in the 100- to 3000-Hz range. In addition, experimentally 
proven hardware fixes, which successfully eliminated the instability, exist for 
each of. the three "correlation" injectors. 

The engine hydraulic stability computer model was run for each of the "corre- 
lation" injectors with both pre-fix and post-fix injector model input. Examina- 
tion and analysis of the model output revealed that the computer model success- 
fully predicted that the fixes applied to each correlation injector would increase 
combustor stability with respect to the instability mode actually observed. 

In Task IV - Sensitivity Analysis, a systematic variation of injector design 
variables was made, using the computer model for the Aerojet and Rocketdyne OME 
technology injectors, to gauge the effects on injector gain. The variables in- 
vestigated included (1) model frequency, (2) orifice resistance, (3) orifice 
inertance, (4) ring groove area, and (5) the presence of dams. For the Rocketdyne 
injector, injector gain was observed to be simply related to orifice resistance, 
orifice inertance, and frequency through a term called the break frequency. 
However, for the Aerojet injector, complications caused by the possibility of 
the pie manifolds acting as Helmholtz resonsators precluded the establishment of 
any simple relationship for the injector gain. Results of these sensitivity 
analyses showed that individual injector design parameters (i.e., orifice char- 
acteristics, ring groove area, etc.) could not be characterized as promoting 
or retarding coupling. Instead, the interaction between all constituent in- 
jector flow passages must be analyzed together. It was concluded, therefore, 
that the development of generalized design criteria is premature to its use in 
place of the parametric application (to a specific injector) of the engine 
hydraulic stability model. The engine hydraulic stability computer model, how- 
ever, was concluded to be quite capable of predicting the relative effects of 
proposed fixes and/or evaluating design options of a particular injector concept. 
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The program was concluded with the conversion of the engineering model from 
Rocketdyne’s IBM 370 computer to the NASA/JSC Univac 1110, EXEC-8 computer sys- 
tem and successful operation of the engineering model at NASA/ JSC. Model docu- 
mentation in the form of the present final report and a computer manual consti- 
tuted the end products of the contract. 

The work performed within all of the foregoing tasks is presented in this 
document. The presentation of the subject matter is organized as a Lask-by- 
task description rather than a detailed discussion of the computer program. 

The latter is extensively described in a separate companion document entitled 
"Engine Hydraulic Stability Model Computer Manual" (Ref. 2). 
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SECTION II 


REVIEW OF EXISTING MODELS 

The development of the fundamental equations describing fluid flow in the in- 
jector is detailed in Appendix A. Depending on the further simplifying assump- 
tions made, these equations can be solved using any of several techniques . 

These techniques form the candidate solution methods that were evaluated in 
Task I. A discussion of each follows. 


COMPILATION OF DATA ON EXISTING MODELS 

The injector modeling techniques which have been employed in existing models 
are described in the following paragraphs . The basic techniques for modeling 
multidimensional wave motion in the constituent portions of the injector are 
(1) lumped parameter approach, (2) continuous parameter approach, and (3) gen- 
eralized wave equations. These techniques were reviewed. The extent of the use 
of these techniques, along with simplifying approximations employed, in existing 
injector models were also examined. Many injector models were developed for a 
"one-shot" analysis of a particular injector. 

Lumped Parameter Technique 

The lumped parameter description consists of dividing any given hydraulic system 
into a number of pressure nodes with flows calculated between those nodes. The 
allowable maximum length between each pressure node is determined by the acoustic 
velocity of the fluid, and the acoustic frequency of interest. This length can 
be explicitly given by the relation 


l itr */« 

where 

L = length between pressure nodes 

a = acoustic velocity of the hydraulic fluid 

v = acoustic frequency 

X = acoustic wavelength 


( 1 ) 
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Equation ( 1 ) states that there must be at least eight pressure nodes per 

acoustic wavelength. 

Figure 5, shows how a given line segment of flow would be broker, into 
lumps using the lumped parameter description. Each line segment consists 
of a series of pressure nodes, P^, separated by a length L. Between these 
pressure nodes are mass flows, w^, which also must be included in the lumped 
description. 


Flow 


r~ l — 1 



Line Segment of Flow 
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Lumped Parameter Description of Line Segment of Flow 


Figure 5 . Lumped Parameter Technique 


The mathematical development of the lumped parameter technique begins from the 
viscous one-dimensional equations for perturbed flow as shown in Appendix A. 
Taking the finite difference of the spatial derivative in Eq. (A-25) y in 
Appendix A, the change in pressure of each pressure node in Fig. 5 can be 
expressed as 
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dP. 
x 

dt 


- w. - w 


i+1 


(2a) 


such that 



(2b) 


O' _ 

w = pAv 

x 

where 

C is the capacitance of the pressure node 
Z is the perturbed mass flowrate 

A is the cross-sectional area of the line segment 
v^ is the perturbed axial fluid velocity 
V is the volume of the pressure node (L times A) 

p is the time averaged fluid density (which is constant) 

P is the perturbed pressure 

t is the time 


(2c) 


and g is a units conversion factor between mass and force 

C , 2 
when Engineering units are used (386 in.-lbm/lbf-sec ) 

Finally, taking the finite difference of the spatial derivative in Eq. (A-26) in 
Appendix A, the perturbed mass flowrate between each pressure node in Fig 5 is 
determined by 


dw. 


/*J 

- - p. - 


dt 1-1 1 


*L*i 


such that 


I = 


A P>, 


*L 


2 Up 


w 


(3a) 


(3b) 


(3c) 
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where 


I I denotes absolute value 

I is the inertance of the fluid element 

is the linearized flow resistance based on 
time-averaged steady-state values 

AP is the time-averaged pressure drop across a 
pressure node of length L 

w is the time-averaged mass flowrate 

Equations (2) and (3) above are the governing equations for perturbed 
steady hydraulic systems, represented by the lumped parameter technique. They 
form a system of linear ordinary differential equations which is often solved on 
an analog computer in real time. 

The lumped parameter technique can, however, be used in either the time domain 
or the frequency domain.* A frequency domain solution is indeed most appro- 
priate when an oscillatory pressure at the injector face is used as input to 
the injector model. In this instance, a Laplace transform of Eqs. (2a) and 
(3a) are taken with the result being 

(4) 

(5) 



* In the time domain, the ordinary differential equations are solved by inte- 
gration as a function of time. In the frequency domain, a Laplace transforma- 
tion is made which converts the differential operator to a Laplace operator. 
This is then replaced by j to and solutions are obtained as a function of fre- 
quency. Frequency response gains and phases are the direct output of the 
frequency domain solution while the output form in the time domain solution 
are individual variables as a function of time. Therefore, to get gains, the 
amplitude and phase of each variable with respect to the input need to be 
determined from the time transients . 



Now replacing "S" by "Jw", where a> is the frequency and j = /-l, a frequency 
response solution (Ref. 3) for the perturbed pressure and mass flowrate at any 
node can be determined by solving the matrix of equations given by Eqs. (4) and 
(5) above. 

The lumped parameter description has previously been used with excellent results 
at Rocketdyne in analyzing several complex injector systems. These include 
the XRL booster, XRL sustainer, and the LE-3 injectors (Refs. 4 and 5). For 
these cases, a specific injector and a specific mode and frequency were ana- 
lyzed. When analyzing a system at only one frequency, computer cost is gen- 
erally small. Computer cost for a given number of equations is directly pro- 
portional to the number of frequencies to be analyzed. 

Continuous Parameter Techniques 

The continuous parameter approach is a method of solution which allows gradients 
of pressure and velocity to exist within the finite segments of fluid flow. 

Thus, the conservation equations describing this situation are partial dif- 
ferential equations, where pressure and velocity are functions of time and 
position. This contrasts with the lumped parameter technique which is charac- 
terized by ordinary differential equations where only time is the independent 
variables, Eqs. (2a) and (3a). 

The continuous parameter technique has had previous usage at Rocketdyne es- 
pecially on the OME f j ed system coupled stability investigation (Ref. 1) > 
and has been shown to have a distinct potential gain in computer core storage 
and accuracy over the lumped parameter technique. This is due to the fact 
that by allowing for gradients to exist in a flow segment less nodes are 
needed to describe the overall flow. This means less equations which implies 
less core storage and less computing time. 

As shown in Appendix A, the development of the continuous parameter technique 
begins from the one-dimensional wave equation for inviscid perturbed flow. 

This representation is found from Eqs. (A-23) and (A-24) of Appendix A which 
give, 
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( 6 ) 


and 


where 
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x = axial direction 


(7) 


The technique used in the OME feed-system coupled stability investigation 
(Ref. 1) to solve the above system of equations, employed a D’Alembert solu- 
tion (Ref. 6). The general form of the solution which satisfied Eq. (6) is 


P = F (t + x/a) + F 2 (t - x/a) 


( 8 ) 


where F^ and F ^ are arbitrary functions. Now, Eq . (7) and (8) can be 
shown to yield 



1__ 

pa 


- ¥ 1 (t + x/a) + F 2 (t - 



(9) 


Figure 8 shows how a long linear flow passage would be broken up into seg- 
ments using a continuous parameter representation. Such segments may re- 
present, for example, passages of various cross sections or passages connecting 
flow junctions. 



Figure 6. Continuous Parameter Technique 
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Using the nodal representation of Fig, 6 above, Eqs. (8) and (9) 
be combined to eliminate the functions, F^ and F^ and show that 





( t— T i ) 



( 10 ) 


and 
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1 A S„ 1 


J ( t~Tf) 
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P . . , t - — w, , . 
i+1 Ag i+1 
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( 11 ) 


where 


= L/a the signal propagation time between modes 


The subscripts, (t ~ x^) and (t) denote the time values at which the terms 
in brackets are to be evaluated. 


Equations (10) and (11) are for inviscid one-dimensional flow. However, 
a viscous one-dimensional flow representation can be included by putting in 
the viscous momentum head loss terms at each node. This is done by simply 

AJ m 

replacing by . Now, taking the Laplace Transform of 

Eqs. (10) and (11), and again replacing "S" by "joj", a frequency response 
solution to the matrix of equations can be obtained as in the lumped parameter 
techniques. 

Multi-dimensional Wave Technique 

This technique allows for a three-dimensional variation of axial perturbed 
flow. Using Eqs,(A-23) and (A-24) of Appendix A will give for the cylindrical 
coordinates of a pipe 


and 
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V'her e 

r = radial direction 
6 = tangential direction 

z = axial diration 

Other cavity geometries may of course be chosen to describe parts of the 
feed system, however, this representation will allow the complexities of 
this method to be seen. 

The simplest solution to the wave equation is obtained when the following 
uniform boundary conditions are given: 


P (r = 0) = finite 


3P 

9r 


= 0 


r=r 


w 


where r = radius of pipe 
w 


n > 

3P 

30 


= 0 


6=0 


(14a) 


(14b) 


(14c) 


and 


P (0=0) = continuous 


(14d) 


Now taking the Laplace Transform of Eqs. (12) and (13) , replacing "S" with 
"jw", and separating variables in Eq. (12) yields 


and 


? - I I 

m=0 k 
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00 
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(A cos k z + B sin k z) cos(m0) 
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i 


where 


C - <^> 2 - (H 


m = 0, 1, 2, 3, ... . 

J = Bessel Function of the first kind, order m 
m 


and the eigenvalues of a are given by the transcendental equation 


d 1~J (a ) ”1 
L m mn J 


dr 


mn 


r=r 


= 0. 


w 


The above constants A and B must still be determined from the boundary 
conditions for the perturbed pressures at the pipe inlet and exit. This 
is a very complex, solution and still no account for viscous momentum head 
losses has been included. 


When non-uniform boundary conditions are given, the solution becomes even 
more complex and requires the use of Green's functions. Some analyses 
using Green's functions were used at Rocketdyne to investigate the potential 
for feed system coupling in the Rocketdyne OME (Ref. 7 ). In that case, 
for frequencies of 2300 to 2800 Hz, simple resonance phenomena were sought 
and only frequencies corresponding to P = 0 at the injector face were 
calculated. 


The complexity of this approach is indeed such that it is way beyond the 
scope of this project. 

EVALUATION OF EXISTING TECHNIQUES 

The configurations of three injectors which have experienced stability 
problems possibly resulting from hydraulic coupling with the combustion 
process were documented in detail. This was done in order to ascertain 
the difficulty or complexity involved in the application of the various 
modeling techniques to describe multi- dimensional wave motion in a "typical" 
injector. Based upon the test evaluation range of variables such as thrust 
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(25 to 5000 lb^) , pressure (100 to 1000 psia) and frequency (100 to 3000 Hz) 
as enumerated in the Statement of Work (see Ref. 8), the maximum injector 
diameter of interest was calculated as shown below. 


D ch 


/ 4(F) (CR) ^ / 4 (5x10^) (2^5) 

J * P c C F J ir(10 3 ) (1.6) 


10 in. 


This 10-inch diameter corresponds to the diameter of a cylindrical chamber hav- 
ing a first tangential acoustic frequency of approximately 3000 Hz. 


Typical Injector Design 

Rocketdyne Lance XRL Booster Injector . Manifold and face pattern details of 
the Rocketdyne Lance XRL booster injector are shown in Figs. 7 through 9. This 
injector comprises an annular area around the central sustainer engine. The 
outer diameter of the annular booster injector is 13.2 inches while the inner 
diameter is 6.65 inches. An unlike doublet orifice pattern is utilized (Fig. 7). 
Orifice diameters are 0.0515 inch for the fuel and 0.073 inch for the oxidizer. 
The XRL injector contains three oxidizer ring grooves and two fuel ring grooves. 
EAch ring groove is fed in four locations from a supply manifold whose shape 
is shown in Figs. 8 and 9. 

Rocketdyne Space Shuttle OME Technology Injector . Details of the Rocketdyne OME 
technology injector are shown in Figs. 10 through 14. The like-doublet orifice 
pattern is shown in Figs. 10 and 11. Orifice diameters are 0.028 inch to 
0.033 inch for the fuel and 0.032 inch to 0.038 inch for the oxidizer. The 
injection orifices are fed by a total of 10 ring grooves (5 ox and 5 fuel) 
behind the 8.2-inch-diameter injector face (see Figs. 10 and 12). The oxidizer 
ring grooves are fed through slanted feeder passages from a central oxidizer 
manifold as shown in Figs. 13 and 14. The fuel ring grooves are fed through 
slanted feeder passages from an annular fuel manifold (see Figs. 13 and 14). 

Aerojet Space Shuttle OME Technology Injector . Details of the Aerojet OME tech- 
nology injector (Ref. 9) are shown in Figs. 15 and 16. The hydraulic diameters 
of the square injection orifices are 0.020 inch for the fuel and 0.024 inch for 
the oxidizer. An 867-element X-doublet platelet pattern on the 8.2-inch-diameter 
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Figure 8. Cross-section of XRL Injector Housing 

Showing Booster Oxidizer and Fuel Passages 
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Section A-A 
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Figure 11. Typical Rocketdyne OME Technology Injector Element 



Typical Fuel Feeder Entrance 
for Fuel Ring Groove 


i 


Acoustic 



tc Ox Ring Groove 


Figure 12. Injector-End View of Rocketdyne OME 
Technology Injector Housing Showing 
Ring Grooves and Feeder Entrances 
to Ring Grooves 
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Fuel Ring Groove (8 Rings) 


Ox Ring Groove (7 Rings) 


Figure 16. Bottom View of Aerojet OMF. 
Technology Injector 
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injector face is fed by a total of 15 ring grooves (8 fuel and 7 ox) as shown in 
Fig. 16. An outer ring circumferential fuel manifold feeds an inner ring fuel 
manifold through 24 drilled holes. Fuel flows from the inner fuel ring through 
three downcomer slots to three pie-shaped manifolds. Downcomer slots from the 
pie manifolds feed each fuel ring. The oxidizer is fed through a central mani- 
fold to three pie-shaped manifolds. As in the case of the fuel, downcomer slots 
from the ox pie manifold feed each ox ring. 

The detail injector drawings show that injectors are, in general, comprised 
of constituent modules which may include (1) inlets, (2) domes, (3) torus 
or ring manifolds, (4) downcomers (holes or slots), (5) pie manifolds, 

(6) ring grooves, and (7) orifices. A myriad variety of geometric varia- 
tions can exit for all of the constituent modules. Domes, for instance, 
may have either single or multiple inlets and often consist of a complex 
geometric shape with usually multiple outlets. A torus may also have either 
single or multiple inlets and outlets. Its geometry may be uniform but is 
often variable so as to result in constant velocity flow throughout the 
torus. Downcomers may consist of drilled cylindrical holes or slots. They 
may be directed radially, axially, or have both radial and axial components. 

Ring grooves may have either single or multiple inlets and always have mul- 
tiple outlets (orifices) . Their geometry may be either constant or variable 
(if constant velocity is desired in the ring groove) . In addition, dams may 
be located in some or all of the ring grooves at particular angular locations. 
The number and size of the injection orifices can vary widely. The orifices 
themselves can compromise many distinctive types of injection elements. 

It was of the utmost importance that the injector modeling technique selected 
for use in the model development task be sufficiently simple and flexible so 
that an injector model could be developed which is both general and yet capable 
of analyzing extremely complex injector geometries. 

Advantages/Disadvantages of Various Techniques 

Three injector modeling techniques (lumped parameter, continuous parameter, 
and multi-dimensional wave solution using Green's functions) were evaluated 
for possible utilization in the generalized injector model to be developed 


in this program. While the three techniques utilize different approximations 
and methods of solution, the basic governing equations for each of these 
techniques is identical. This is shown in detail in Appendix A. 

The advantages and disadvantages of the afore-mentioned injector modeling 
techniques which were initially apparent are summarized below. 

Lumped Parameter Technique 

• Advantages 

• State-of-the-art 

• Simplicity 

• Extensive previous usage 

• Disadvantages 

• Core storage 

• Computer time 

• Potential limitation on injectors which can be described adequately 
with core storage limits 

• Potential accuracy due to possible description limitation 

Continuous Parameter Technique 

• Advantages 

• Previous usage 

• Potential gain in core storage and computer time accuracy 

• Disadvantages 

• State-of-the-art does not exist for application to line segment which 
includes mass gains or losses along its length (i.e., a ring groove 
segment feeding injection orifices). 

Multi-Dimensional Wave Solution (Green's Function) 

• Advantages 

• Inherent high frequency capability 

• Computer core storage 

• Disadvantages 

• Difficulty in handling complex boundaries 

• Difficulty in handling intermediate flow junctions 

• Mathematical instability of solution 
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Irrespective of the technique selected, the most appropriate output from the 
injector model should be gain and phase as a function of frequency relating 
injector flowrate as a function of the chamber pressure perturbation. Flow 
and pressure distribution throughout the injector is thus determined. In 
order to obtain the model output in this form, the governing injector model 
equations are linearized and subsequently arranged in matrix form. The coef- 
ficient matrix and input matrix serve as input data to the frequency response 
program, which then yields the required injector frequency response (see sche- 
matic representation below) . 



This approach was selected after careful consideration for use in the OME feed- 
system model (Ref. 1).* 

Use of the frequency response program is cost effective since time transients 
are not included. Thus, a "steady-state" oscillatory determination of pres- 
sure and flowrate is obtained. The frequency response method is thus generally 
preferred over more lengthy (and costly) solutions in the time domain (Ref. 10)* 

Lumped Parameter Technique 

The ability of the lumped parameter technique to adequately describe the com- 
plex injectors of interest within core storage and computation time limita- 
tions was deemed an item of urgent concern. 

As stated earlier, the advantages of the lumped parameter description are 
that the technique is simple, flexible, and has extensive previous usage in 
analyzing the complex injector flow patterns. The disadvantages are the 
tendency to require larger core storage and computer time in order to have 
the injector adequately described. A promising technique was investigated 


* Further, the statement of work called for an injector model whose structure and 
format is compatible with this existing generalized OME feed-system model. 
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to eliminate these disadvantages. An explanation of how the contemplated 
method would be used on a typical ring groove (Fig. 17) is discussed below. 


IN 



PC(1) 


PC(2) 


PC(3) 


PC(N) 


Figure 17. Typical Ring Groove for Evaluating Method 


As shown in Fig, 17 , for a ring groove with N pressure nodes, N pressure 
equations and 3N flow equations are required. Therefore, for 12 pressure 
nodes, 48 equations are required. Since each variable and the input re- 
quire a real and imaginary term, the matrix set up in the frequency response 
program would have to be a 48 by 98 matrix. For a complex injector with 
several ring grooves, the size of the matrix computer time for inversion 
would become very large. The method under evaluation involves the elimina- 
tion of all the flows from the set of equations by substitution. As an 

a* 

example, the equations for P(2) are: 

SP (2) = ^ — |_*(2) + w(N+2) - w(N+3) - *(2N+2)J 

Vg c 

where 

S = d/dt 
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1(2) S w(2) = P IN - P(2) - R(2) w(2) 

I(N+2) S w(N+2) = P(l) - P(2) - R(N+2) w(N+2) 

I(N+3) S w(N+3) = P(2) - r(3) - R(N+3) w(N+3) 

I(2N+2) S fr(2N+2) = P(2) - PC(2) - R(2N+2) w(2N+2) 


where I 


L 

Ag. 


Each flow equation can be solved in the form 


w(2) = 


P IN - P(2 > 
R(2) + 1(2) S 


The flow equations can then be substituted into the pressure equation to give 


S P(2) 



P — P ( 2 } A* 

^IN P(l) - P(2) 

R(2) + 1(2) S R(N+2) + I(N+2) S 


P (2) - P(3) P (2) - PC(2) 

R(N+3) + I(N+3) S “ R(2N+2) + I(2N+2) S 


By multiplying each term on the right hand side by R - IS, and substituting 
j to for S, the equations will contain only pressure variables and be in the 
correct form for the frequency response technique. Thus, a system with 12 
pressure nodes would be only 12 equations and the required matrix size 
would only be 12 by 26. This means that four times as many pressure nodes 
can be selected and require the same amount of computer time as the fre- 
quency response would require before substitution. It has been determined 
(by the comparison of results using single precision and double precision 
that single precision gives answers almost identical to those obtained using 
double precision, therefore, double precision is not required. 

The above solution technique significantly reduces the core storage require- 
ment for the lumped parameter technique. Thus, more pressure nodes can be 
selected (and the injector consequently more accurately described) without 
increasing computation time. While this improvement does much to eliminate 
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the lumped parameter technique disadvantages of a larger core storage 
requirement and increased computer time, it was felt absolutely necessary 
to demonstrate that a "typical" injector could indeed accurately be des- 
cribed by about 100 pressure nodes,* 

Of the injector configurations presented in Figs. 7 through 16 , the 
Aerojet OME technology injector is believed to be the one that will require 
the largest number of pressure nodes to accurately define the system be- 
cause of the large number of annular manifolds and ring grooves. A schema- 
tic of the lumped parameter distribution for the fuel side of the Aerojet 
OME injector is shown in Fig. IS . The means in which the constituent 
modules are divided into pressure lumps (or nodes) is shown in Fig. 19.** 
Each box in Fig. 18 represents a pressure node and each line between pres- 
sure nodes represents a flowrate. This lumped parameter model thus repre- 
sents all the constituent injector modules including the outer and inner 
manifold rings, the pie manifolds, the ring grooves, all the flow passages 
between the manifolds and ring grooves, and chamber pressures. This lumped 
parameter distribution has 99 individual pressure nodes and 223 individual 
flowrates. A frequency response model of this size will easily fit on the 
IBM or Univac computer. 

The determination of this system of pressures and flowrates was based on 
the use of eight elements (nodes) per wavelength. At a frequency of 3000 
Hz, each pressure node thus has a length equal to or less than approximately 
2.25 inches. This criteria, while definitely adequate for accuracy (as 
shown later), may not be necessary. A discussion of the analysis of the 
results of using different numbers of lumps on a single ring groove is 
presented later in this section of the report. 

In order to run the generalized injector computer model, a number of inputs 
need to be specified. The lumped parameter distribution which accurately 

*It is believed that the core storage requirement on the Univac computer 
should allow at least 140 pressure nodes . 

**The inner ring manifold is displaced from its actual location in Fig. 19 
for clarity of understanding. 
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Outer Manifold Ring (15 lumps) 

Inner Manifold Ring (15 lumps) 

•Ring Groove No. 8 (12 lumps) 

Ring Grooves No. 6, 7 (9 lumps apiece) 

Ring Grooves No. 4, 5 (6 lumps apiece) 



Plane of ring 
groove dams (typical) 


Pie Manifold 
(6 lumps) 


Ring Grooves No. 1, 2, 
(3 lumps apiece) 


Figure 19. Distribution of Lumps for Fuel Side of 
Aerojet OME Technology Injector 
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describes the system (like Fig. 18 ) must be selected. Several options 
are available for generating the system of equations required to util- 
ize the frequency response program. The system of equations consists of 
a pressure equation for each node and a flow equation for each flow between 
pressure nodes. Each linearized equation for pressure can be written (using 
Eq. (2a) as its source) as: 

r\> 'i w n> 

P = (E w. - Z w ) 

S Vg xn out 


P 


W. 

xn 


out 

S 


V 



= peak-to-peak oscillatory pressure, psia, 
= oscillatory flow into element - Ib/sec 
= oscillatory flow out of element - lb/sec 
= Laplace operator 

3 

= fluid volume in pressure node - in. 

= 386 lb in. /lb, sec^ 

m f 

= capacitance of fluid element 


Each linearized equation for flow can be written (using Eq. (3a) as its source) 
w - As c (P - P - M) 

ST u B 

where 

w = peak-to-peak oscillatory flow - lb/sec 
L = length of flow element - in. 

S = Laplace operator 

2 

A = flow cross-sectional area - in. 

g = 386 lb in. /lb- sec^ 

a c m f 


as : 
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rj 

P 


U 


D 

R 


L/Ag„ 


peak-to-peak oscillatory upstream pressure, psia 

peak-to-peak oscillatory downstream pressure, psia 

_ _ 2 

linearized flow resistance = 2 AP/w, sec/in. 

2 2 

fluid inertance, lb,. sec /lb in. 

f m 


For these equations, the inputs consist of effective fluid acoustic velocity, 
volume of each pressure, node, fluid inertance and linearized resistance for 
each flow equation. In addition, the logic for coupling all the flow and 
pressure equations together with each chamber or upstream pressure input is 
required. For some typical injectors which fit a specific lumped parameter 
configuration, the logic for coupling the flows and pressures together may 
be contained in the program. For any configuration that does not fit a 
typical set of equations, the coupling must be included as part of the input. 
The flows entering and leaving each pressure node must be directionally speci- 
fied as well as which flows terminate into each of the different input chamber 
pressures. As an example, take one of the ring groove pressure nodes that 
might have a flow distribution as follows: 
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For this pressure node, the inputs would be the number of flows entering 
or leaving the pressure node (four in this case) , and the subscript of 
each flow (+k and +j for flows entering and and -m for flows leaving) . 

In addition, the subscript of each flow terminating in each chamber pres- 
sure input must be specified. The amplitude and phase (0 or 180 degrees) 
of all the input pressures must be input referenced to one specific loca- 
tion. From this data, in addition to the inertias, resistances, and capa- 
citances, the computer program could set up the equations, do the matrix 
manipulation, and print the output results. 

The output of the model is gain and phase of each dependent variable (all 
the individual pressures and flowrates) with respect to the referenced inlet 
pressure. This provides the user with the complete flow and pressure distri- 
bution throughout the Injector. 

The lumped parameter technique could be solved in either the time or the 
frequency domain. The advantage of the time domain is that nonlinearities 
can be included. However, there are several disadvantages of the time 
domain. These include: 

1. Computer time required 

2. Insuring convergence of solution technique 

3. Determination of when results have reached 
a constant oscillatory amplitude. 

To calculate the injector operation in the time domain, all the equations 
are written, and the input amplitude and frequency are specified. Then the 
injector computer model is allowed to run for as many cycles as required to 
allow all the variables to reach a repeatable (constant) oscillatory ampli- 
tude and phase relationship. For a complex model, very small integration 
time intervals are generally required to keep the system of equations digit- 
ally stable. The net result is usually a large expenditure for computer 
time. This large cost is one of the main reasons that the frequency response 
technique was originally developed. Another reason is that the frequency 
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response solves the equations directly with no iteration or step type cal- 
culations required. Therefore, for linear systems analysis, the frequency 
response technique should be used. 

The lumped parameter technique, using a frequency response method of solution, 
was used to analyze a typical ring groove as described in Fig. 17 . This 
method -of solution (described earlier) was verified by comparison to a 
standard frequency response method* and then used to evaluate the effect of 
various number of pressure nodes in the ring groove. The ring groove analyzed 
had a total length of about 30 inches. Based on the criteria of eight elements 
per wavelength, this system could be accurately described by about 14 pressure 
nodes. The input chamber pressure profile corresponded to a chamber first 
tangential mode. Table I shows a comparison of the response of the ring groove 


TABLE I 

GAIN OF RING GROOVE PRESSURE TO CHAMBER REFERENCE PRESSURE 


Frequency 

0 Hz 

2500 Hz 

8000 Hz 

Number 

of Lumps N. 




6 

0.0869 

1.382 

0.1977 

12 

0.08304 

1.35 

0.1984 

24 

0.08209 

1.341 

0.1986 

48 

0.08185 

1.3394 

0.1986 

96 

0.08183 

1.3392 

0.1986 


*The standard frequency response method solves the full set of pressures and 
flowrate equations without doing any manipulation to eliminate variables. 
This method has extensive use and verification at Rocketdyne. The method 
of eliminating variables prior to the matrix inversion reduces the computer 
time and storage requirements. The test case run was compared to the stand- 
ard frequency response to verify that the accuracy of the matrix solution 
process w 7 as not affected by the alteration of the equations. 
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pressure directly above the reference pressure for various numbers of pressure 
nodes selected. (Naturally, at zero frequency it would not be possible to have 
a first tangential chamber mode, but results are included for comparison pur- 
poses.) These results show that even using only six lumps gives results that 
are within 6.2% of the actual frequency response* at zero Hz, 3.1% at 2500 Hz, 
and only 0.5% at 8000 Hz. Two significant conclusions may be drawn from these 
results. First, the use of the crtieria of eight lumps per wavelength will 
definitely provide good results since any frequency response within 6% is gen- 
erally satisfactory. Second, the technique gives accurate results for fre- 
quencies higher than 3000 Hz as long as the first tangential mode is being con- 
sidered. This is probably due to (1) a reduced amount of response in the feed 
system at higher frequencies,, or (2) a reduction of the distance into the feed 
system into which the input can propagate at higher frequency. It may be con- 
cluded that this lumped parameter frequency response technique is not limited 
to 3000 Hz, but may be accurately applied to higher frequencies as well. 

Continuous Parameter and Multi-Dimensional Wave Solution Techniques 

The continuous parameter and multi-dimensional wave solution (Green's function) 
techniques initially appeared attractive because the division of the injector 
into a nodal network (as is used in the lumped parameter technique) is not 
required. (The injector still may be divided into flow segments but the fur- 
ther division of these flow segments into nodes or lumps is unnecessary.) This 
possesses a potential advantage in terms of core storage and computation time. 
The extremely complex and nonuniform boundaries existing between the constituent 
modules of a "typical" injector present, perhaps, an even greater obstacle. 

Examination of the continuous parameter technique led to the conclusion that 
the state of the art at this time does not exist for the application of this 
technique to line segments which allow' for mass gains and/or losses along 
its length (i.e., a ring groove segment feeding injection orifices). An attempt 
was made to modify the continuous parameter technique such that it would be 
capable of describing flow segments having mass loss or gain. 


* The actual frequency response is assumed to be that obtained using 96 lumps. 
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The conservative equations of mass and momentum were written for cylindrical 
control volume (considered to represent one line segment in a ring groove). 
While partial differential equations were developed and appropriate initial 
and boundary conditions established for describing flow segments having mass 
gain and/or loss with the continuous parameter technique, the resulting system 
of equations were not adapted to a frequency response type of solution. The 
progress made, therefore, was not sufficient to bring the continuous parameter 
technique to the point where a "side-by-side" computer comparison with the 
lumped parameter model could be made. 

For feed lines to the ring grooves, waterhammer (continuous parameter) equa- 
tions could be used, if this would be of value, even if the ring grooves were 
described by a finite difference network. The current Rock.etdyne frequenpy 
response program (used in the hydrodynamics section of the OME Feed System 
Coupled Stability Investigation, NAS9-14315 (Ref. 1)) routinely handles com- 
binations of waterhammer equations and finite difference equations. However, 
for "typical" injector designs such feed lines are usually short enough (less 
than 1/8 wavelength) to not require this type of presentation. Indeed, the 
combination of the continuous parameter and lumped parameter techniques to 
describe a single injector may increase the amount the core storage required. 

The possible application of the Green's function technique to this program 
was discussed with Dr. Carl Oberg (who was responsible for all previous ap- 
plication of this technique at Rocketdyne) . It was agreed that, based on 
current knowledge, the application of this technique to the injector model 
possesses little likelihood of success within a reasonable amount of time and 
dollars . 

Finally, the three candidate modeling techniques were evaluated against a 
number of criteria as shown in Table II. This evaluation was done (for each 
technique) in respect to a "typical" injector configuration whose character- 
istics have been discussed earlier in this report. 



TABLE II. EVALUATION OF CANDIDATE TECHNIQUES 


" ' — -^TECHNIQUE 
CRITERIA ' 

LUMPED 

PARAMETER 

CONTINUOUS 

PARAMETER 

GREEN'S 

FUNCTIONS 

State-of-art 

.. . ... . 

Within state-of-art 

Not within state- 
of-art 

Not within state- 
of-art 

Previous Usage 

Extensive 

Limited 

None 

Accuracy 

Sufficient 

Sufficient 

Sufficient 

Complexity 

Complex 

More complex 

Most complex 

Computer storage 

Acceptable storage 
required 

Acceptable storage 
required 

Least storage 
required 

Computation, time 

1 • 1 ; 

Acceptable (1 min.) 

Acceptable 

Longer 

Type of input 

Physical geometric 
parameters 

Physical geometric 
parameters 

Physical geo- 
metric parameters 

Ease of input 

Laborious 

Laborious 

Impractical 

Type /us ability of output 

Numeric/good 

Numeric /good 

Numeric/good 

Number & obtainability 

Gain and phase/ 

' Gain and phase/ 

Gain and phase/ 

of characterization 
parameters 

adequate 

adequate 

adequate 

Numerical stability 
of solutions 

Not a factor 

Not a factor 

Definite problem 

Applicability to 
digital solutions 

Good 

Good 

Good 

Limitations 

Detail of input 
specification 

Detail of input 
specification 

Detail of input 
specification 























































SELECTION OF MODEL TECHNIQUE 


In evaluating how the various candidate techniques could be applied to 
"typical" injector configurations (see Table II), it becomes apparent that 
the flexibility and versatility of the lumped parameter technique offers a 
great advantage. 

This technique is extremely well suited to adequately describing even the 
most complicated injector configurations. Injectors having a lack of sym- 
metry do not present a problem. The presence of ring groove or manifold 
dams can easily be handled by assigning a very large resistance to the 
particular flow branch thereby effectively eliminating its presence in the 
injector system (Ref. 1), Disadvantages previously believed associated with 
this technique (i.e., core storage, computer time, adequate injector des- 
cription, high frequency capability) are believed to have been resolved. 

In contrast, the other candidate techniques have been shown to require 
considerable advancement of the state-of-the-art before the feasibility 
of utilizing these techniques to model a "typical" injector could be 
satisfactorily determined. 

The selection of the lumped parameter technique for utilization in the 
model development task thus appears to be the only reasonable, cost-con- 
scious choice. 
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SECTION III 


MODEL DEVELOPMENT 

A digital computer model of the injector manifolding which describes the 
coupling of the injector hydraulics with the combustion process was formulated. 
This model was designed to be compatible with an overall generalized engine 
system dynamics model, developed by Rocketdyne for NASA/ JSC under contract 
NAS9-14315 (Ref. 1), which includes propellant feed system hydrodynamics, com- 
bustion dynamics, and chamber dynamics. Compatibility is derived from the use 
of the complex matrix solution technique in both models. This technique can 
use either lumped parameter or continuous (distributed) parameter equations. 

The two models could be coupled by properly setting up the complex matrix with 
no changes in the equation types from either model. The generalized engine 
system dynamics model uses a single lump to describe the injector, while the 
engine hydraulic stability model employs a multi- lump (multi-dimensional) in- 
jector description. The engine hydraulic stability model cannot be. coupled 
into the generalized system dynamics model at the present time because the 
latter model does not employ a multi-dimensional description of the combustor 
and combustion process. 

FORMAT AND STRUCTURE 

Generalized criteria which were employed to develop a format and structure for 
the injector model are as follows: 

a. Input and output formats shall be appropriate to allow incorporation 
of the newly developed injector /manifold model in the generalized 
system dynamics model developed under contract NAS9-14315. 

b . The model shall be structured to quantitatively evaluate the influ- 
ence of injector design variables and geometry in terms of resist- 
ance, capacitance, and inertance on the ability to hydraulically 
couple with the combustion process. 

c. The model shall be capable of analyzing a single portion of the 
injector, such as an orifice, propellant channel, downcomer, etc., 
independent of the remainder of the injector as well as be capable 
of analyzing the entire injector. 
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d. As a minimum the model shall accept an input pressure profile 
from either the chamber or the injector inlet. Profiles 
selected are as follows: For the inlet, a uniform pressure at 

any frequency for 100 to 3000 Hz. for the chamber, either a 
uniform profile or an acoustic profile corresponding to a parti- 
cular chamber mode with nodal diameter orientation and frequency 
ranging from 100 to 3000 Hz. 

e. The model output shall include a single-value gain to simplify 
model utilization and interpretation. 

f . The model shall be formulated as simply as possible consistent 
with the quantitative sensitivity and shall require a minimum 
of engineering judgement factors. 

The engine hydraulic stability computer model uses a frequency response 
program to solve a matrix of simultaneous linear equations. Since this 
is the same basic program used in the generalized engine system dynamics 
model, developed under contract NAS9-14315, the input and output formats 
would allow incorporation of this model into the generalized model. 

As an example of showing how the model format and structure is developed, 
consider the fuel side of the Aerojet OME technology injector as described 
in Section II. Each of the pressure nodes, flow paths, and input pressures 
are numbered sequentially as shown in Fig. 20. This system has 99 dependent 
pressure variables (denoted by rectangles), 52 input pressures (circles) and 
223 flowrates (arrows) . 

This feed system description provides the format for the introduction of the 
input parameters and solution technique. The following discussions show how 
the input parameters to the model are determined and how they are fed as in- 
put into the model. 
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MECHANIZATION OF DETAILS 

As discussed previously, an injector is described in terms of which flow seg- 
ments enter and leave each pressure node. For each flow segment, a resistance 
and an inertance are required as model inputs (line lengths and areas are used 
in determining inertance). For each pressure node, volume and acoustic velo- 
city are required as model Inputs (volume and acoustic velocity determine 
capacitance) . The following paragraphs describe the methods used to calculate 
specific values for these model inputs. 


The model input parameters which need to be determined when using the lumped 
parameter technique are capacitances, inertarsces, resistances, and input 
pressures (amplitude and phase) . The criteria for calculating these model 
input parameters are given in the following discussion. 


Capacitance 

The fluid capacitance term is developed in Section II and given by Eq. (2b) 
as 



V i s c 


(2b) 


where 


C. 

1 

V. 

x 


8 C 


a 


is the capacitance of the i pressure 

is the volume of the i^ pressure node 

is a units conversion factor between mass and force 
when engineering units are used (386 in-lbm/lbf-sec^) 

is the. acoustic velocity of the hydraulic fluid 


In calculating the volume of a pressure node, its volume must also include 
half the volume of each flow segment connected to that particular pressure 
node. Thus, the summation of all nodal volumes must equal the total volume 
of the injector. 
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Inertance 


The fluid inertance term is also developed in Section II and given by Eq. (3b) „ 
However, when flow segments of multiple cross-sectional areas are used Eq. (3b) 
must be changed to 


I. 

1 



where I, 


til 

is the inertance of the i flow segment 

is the length of the flow passage within 
til 

the i flow segment 


A. 

3 


cross-sectional area of the j flow passage 
within the 1^ flow segment. 


(17) 


The inertance equation shown above takes into account that each flow segment 
may be composed of line lengths which have different cross-sectional areas. 
The overall length of one flow segment is determined by the actual distance 
between one pressure node center and another where the flow can be entirely 
characterized as one-dimensional. 

Resistance 

The linearized flow resistance is developed in Section II and given by 
Eq. (3c) as 

2IAPJ 


is the linearized flow resistance of the i C ^ 

flow segment based on time-averaged steady-state values 

is the time averaged steady-state pressure drop 
across the i*"* 1 flow segment 


Where 


V 


AP. 
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w is the time-averaged mass flowrate through the 
.th fl 

i flow segment 


If pressure drop and flowrate, data were known for each flow segment, cal- 
culation of the linearized flow resistances, 's becomes an easy task. 
However, when such data is not known the time averaged steady-state values 
of AP^ -and w^ must be calculated analytically. The procedure used for 
analytically finding the pressure drops and flowrates is discussed below. 


(18) 

where 

Kj; is the pipe friction head losses (4f L/D) 

is the minor friction head losses due to entrance, exit, 
elbows, orifices, etc,, within the line flow segment 

v is the time-averaged steady-state fluid velocity within 
the i*"* 1 flow segment 


A P . = 

1 


K f + 


11 

2 

m=l 


K. 


m 


2 

_£JL_ 

2 S„ 


th 

It should be noted here that if the i flow segment consists of various cross- 
sectional areas all head losses must be referenced to the flow through one part 
of the flow segment. For example, consider Fig. 21 below which consists of a 
flow segment with line lengths of cross-section A^, A ? , ... A. respectively, 
and velocities, v^, ... respectively. 
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Some of the K's will be referenced to line 1, others to line 2, and still 

others to line j . Therefore, AP . is written as 
J ’ 1 


AP. = 


i 2 8, 


(K- + EKJ v." r (K, + E K ) v, 
f t 1 t t t 


+...+ (K f + n t ) v 2 

m . 

J 


However, using continuity, it is known that 


w. = (p A. V- ) = (p A_ v„) = (p A. v.) 

i i i JJ i 


Referencing all head losses to lire j by substituting Eq. (20) into 
Eq. (19) yields 


AP. = R. w. 

l ii 


where 


R. = 

l 


(K f + E K t ) 2 + (K f + Z K t } 2 + ' • * + (K f + E K t > 

m 1 A 1 m 2 A 2 

0- A 2 

2p A j g c 


Any injector can be described by its R ''s regardless of total pressure drops 
and flowrates, since R^ is determined by geometry only. Therefore if an 
injectors R^’s can be combined and reduced to one overall injector resistance, 
R^,, then Eq, (21a) shows that total injector pressure drop, AP^, can be deter- 
mined directly from total injector mass flow, w^, by the relation 


- — 2 

AP t = R t w t 
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Usually total injector mass flowrate is known so that once R^, has been found 

from the individual RAs, then the total pressure drop can be determined 

directly. Once AP m has been found from w„, all individual AP.'s and w.'s 
T T ix 

can be determined by separating R back into its individual components 
analogous to electric circuit resistance reductions. These circuit reduc- 
tions for fluid flow have been developed for series and parallel flows and 
are described below. 


First consider i flow segments in series with each other as shown in Fig. 22. 


AP, 


AP, 


w, 


T ' 


o — VV — ° — w 


-O— • • • — o- 


AP. 

A V 

R. 

x 


Figure 22. Series Flow 


It is known that the total pressure loss through all i flow segments 
AP tS , is given by 


AP ts = AP 1 + AP 2 + . . . AP 


(23) 


or 


A? TS = R 1 *1 2 + R 2 *2 + • • • R i 


(24) 


From continuity it is known that 


*1 = w 2 = w. = * Tg 


(25) 


where 


w is the total flow through the series flow segments. 
J. b 


Therefore the total pressure loss across the series flow circuit is found 
from Eqs. (24) and (25) 


As - < R 1 + R 2 + • • • V *TS 2 


(26) 
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Now consider i flow segments in parallel with each other as shown in Fig. 23 




Figure 23. Parallel Flow 


It is known that 


w T p = Wp + w 2 + ... Wp. 


(27) 


Also the pressure drop across each flow segment is equal, thus 


AP xp = AP X = AP, = AP. 


(28) 


where AP^p is the pressure drop across a parallel flow circuit. 

Finally, using Eqs. (21a), (27), and (28), it can be shown that for flow 
segments in parallel, 

2 

- 2 

W TP 

(29) 
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Therefore, if an injector can be reduced to a simple network of flow resistances 
in series and parallel with one another, then Eqs. (21) - (29) can be used to 
analytically solve for the AP^s and ' s . This in turn allows for the solu- 

tion of the ' s which are then inputted into the computer model. 

The fuel side of the Aerojet OME technology Injector, completely diagrammed in 
Fig, 20j was reduced to a simple network of series and parallel flows as illus- 
trated in Fig, 24. This network reduction was obtained by assuming negligible 
pressure drops in flow segments 2-16, 32-46, 53-73, and 122-172 of the OME fuel 
side manifolding. The network resistances used in Fig. 24 were obtained from 
Eq. (21b) using typical criteria for the K^'s and K^'s as found in the SAE 
Aero-space Applied Thermodynamic Manual (Ref. 11). Using the exact geometry 
of the injector (as reported by Aerojet in Ref. 9), a total injector pressure 
drop of 54.8 psid was analytically calculated for a fuel mass flowrate through 
the injector of 7.19 lbm/sec. This compares very well to the 57 psi pressure 
drop reported by Aerojet. 


Input Pressure 

For any feed system coupled instability above a frequency of 1000 hz, the com- 
bustion chamber response comes either from a combustion chamber resonance or 
from response of pressure in the cup of a recessed post injector. To analyze 
the case of a chamber resonant mode with the injector stability model, the 
relative amplitude and phase of the chamber pressure oscillations across the 
injector face must be input to the model. 

An analytical formulation describing the perturbed flow in the combustion 
chamber was developed for the types of chamber instability modes commonly 
experienced during hot firing. This was especially needed in knowing how 
to input a standing tangential combustion chamber mode as opposed to a spin- 
ning tangential mode. 
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Figure 24 . 


Simplified Resistance Network for Fuel Side 
of Aerojet OME Technology Injector 



This analytical formulation, in the time domain, is shown in Appendix B. 

However, since the engine hydraulic stability computer model employs a 

frequency response solution, pressure inputs must be expressed in terms 

of gains and phases. For tangential pressure modes it is convenient to 

reference all pressures to the pressure at the combustion chamber location 

of r = r , and 0-0 (where r, r^, and 0 are defined as in Appendix B) . 

The standard procedure in Frequency Response Theory (Ref. 12) is to define a 

new variable 0 such that 
m,n 


6. 


m,n 


rf? } 

ot- 1 m,nj 


^{ Pm ’ n ref} 

(30) 


where P 


m ’ n | ref 

the Laplace transform defined as 


is the perturbed pressure at r = r , and 0 = 0, and a£. 


is 


2 ( ? »J - / 


-S t 
m,n 


P dt 
m,n 


where t = time 


Thus using Eq, (30) above, the time dependency can be removed from Eqs. 

(B-12) and (B-18) of Appendix B. By replacing "S " with "jto " (where 

m,n xn y n 

j = r=r and n = frequency of oscillation in radians/sec) the frequency 


input for a standing tangential wave is 

/to r \ 

j 

R = E 1 -is / 


cos m 0 


m,n 


( to r 
a 


(31) 


and for a spinning tangential wave 


(o r 
m,n 


m \ a. 


3. 


m,n 


f.. r \ 

‘Mj ; 1 


cos m 0 + j sin m 0 


m \ a 


(32) 
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where in 


0 12 3 


n 1 ) 2, 3- • • « 

a c = sonic speed in chamber 

J = Bessel function of the first 
kind order ra. 

Equations (31) and (32) above then give the needed information for expressing 

t'he chamber pressure input in terms of gains and phases. . This is because 

£ is a complex vector on the S -plane. The magnitude of this vector 
m,n m,n 

is known as the gain, and the angle it makes with the positive real axis is 
known as the phase. Graphically this is shown in Fig. 25 below. 

Imaginary ( j ) axis 



From Eq. (31) and Fig. 25 it can be seen that standing waves have only phase 
shifts of 0 and 180 degrees depending upon whether the expression on the right- 
hand side of Eq. (31) is positive or negative. On the other hand Eq. (32) 
shows that spinning waves will have phase shifts over the entire range of 0 
to 360 degrees. 

Figure B-3 of Appendix B shows the normalized maximum amplitudes for the first 
and second tangential modes and for the first iradial chamber mode as a function 
of the non-dimensionalized radius. To simplify the input to the model (the 
fuel side of the injector has 51 separate chamber pressure inputs) , curve fits 
were made for these three modes. The curve fits used are: 
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IT Mode - standing 

Amp = (1.5056 R + .0641 R 2 - .5697 R 3 ) COS(0) 

IT Mode - spinning 

Amp = (1.5056 R + .0641 R 2 - .5697 R 3 ) 

2T Mode - standing 

Amp = (-.17 R + 3.583 R 2 - 2.413 R 3 ) COS(29) 

2T Mode - spinning 

Amp = (- . 17R + 3.583 R 2 - 2.413 R 3 ) 

1R Mode - 

Amp = 1 + .15R - 4.87 R 2 + 3.315 R 3 

where 9 is the selected injector tangential location and R is the normalized 
radius. The location where 6 = 0 is defined as that location where the phase 
is zero degrees and the amplitude of a tangential input is maximum. The re- 
quired model inputs to input one of these chamber modes are the radius and 
angle for the location of each chamber pressure input, and whether the mode 
is standing or spinning.* The model then calculates the amplitude and phase 
inputs for the selected chamber pressure mode. The model also has the option 
for any general input matrix. . That is, the programmer can specify the ampli- 
tude and phase at each chamber pressure input. 

INPUT AND OUTPUT ROUTINES 

The input and output formats used in the model are compatible with those 
employed in the feed system coupled model developed by Rocketdyne for 
NASA/ JSC under contract NAS9-14315 in the event that incorporation of 
the injector model into the latter model be desired at some later date. 

Input routines were written which allowed the selection of either an imposed 
feed system pressure profile or an imposed chamber pressure profile (either 
uniform or characteristic of a particular chamber acoustic mode) with an 
arbitrarily selected attendant oscillation frequency. Output routines 
consist basically of the determination of flowrate and pressure variations 

* A discussion of standing and spinning waves is found in Appendix B. 
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throughout the injector as a function of frequency for each imposed input pres- 
sure profile. 

To provide a simple way to evaluate any specific injection configuration, out- 
put routines were included which allow printout of various summations of the 
oscillatory injector flow (for a unit pressure perturbation). These summations 
include (1) a summation of all the absolute values of injector flowrate, (2) a 
vector summation of all injector flowrates (attempt to include phase angle), 
and (3) and (4) the previous two cases except that each injector flow is multi- 
plied by the relative amplitude of the chamber pressure that it feeds before 
making the summations. These summations are .given in both lb/sec/psi and per- 
cent chamber pressure. 

The first summation simply adds all the gains for all the injector flows. 

This provides a measure of the total injector flow response with no considera- 
tion as to how the flow could couple with the combustion process. The second 
summation considers the phase relationship between each injector flow and the 
chamber pressure flow that it feeds. If, for example, two flows that were 
180 degrees out of phase both fed chamber pressures that had the same phase, 
it is unlikely that both flows could cause coupling since one or the other 
would have the wrong phase relationship. Therefore, the second flow summation 
was made where the flows were vectorially added . The phase assigned to each 
flow is based on the difference between the phase of the floy and the phase of 
the pressure that it was feeding. 

Another characteristic that can affect the amount of coupling is the phase of 
the combustion chamber response. A flow that is feeding a location in the 
chamber that is a pressure node cannot cause coupling, whereas a flow feeding 
a pressure anti-node location would have maximum coupling potential. There- 
fore, the third and fourth summations were made where the amplitudes of each 
flow was multiplied by the respective chamber pressure amplitudes before 



summing the values. All four summations are printed out so that the user can 
evaluate each of them.* 

In addition to printout of the summation of injector flows (both absolute and 
vector) , the ease of interpreting the output of the engine hydraulic stability 
computer model was also enhanced by an optional graphical display of the ampli- 
tudes of ring groove flows, ring groove pressures, injector flows, and chamber 
pressure inputs around each ring groove. 

MODEL CHECKOUT 

In order to first model the system the computer program must be told how many 
dependent and input pressure nodes there are. Also, the program must be told 
what flows are associated with each pressure node and what pressure nodes are 
input nodes. This tabular computer input is shown in Fig. C-l of Appendix C. 
Flows assumed to enter the node are denoted by a plus sign and each flow assumed 
to leave a node by a minus sign. Each pressure node can have up to 12 flow 
terminations (6 entering flows and 6 leaving flows). If the node has less 
than the maximum number of flow terminations, the additional flows are read in 
as zero. For pressure node 1, for example, flow 2 enters, and flows 3 and 17 
leave . 

The model was run for the fuel side of the Aerojet OME technology injector 
using the resistances, capacitances, and inertias previously calculated, and 
using a spinning first tangential mode input at 2600 Hz (frequency of observed 
instability). The method for determining these resistances, capacitances, 
inertias, and pressure inputs was the same as that outlined in the "Mechaniza- 
tion of Details" portion of this section of this report. 


* Generally, on the model simulations made so far in this program, all four 
summations tend to change in the same direction for any specific injector 
modification so that no one summation has to be considered the correct one 
to use. If one did have to be picked, the best one would probably be sum- 
mation number four that is the vector sum times the chamber amplitudes. 
However, in a case where a modification had different effects on the dif- 
ferent summations, any conclusions made should be considered questionable. 
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Appendix C gives a tabulation of the model input and output data for the fuel 
side of the Aerojet injector. A list of the input data deck used for this 
run is shown in Fig. C-l. Fig. C-2 shows the amplitude and phase inputs cal- 
culated in the program. For each flow, the printout shows the upstream and 
downstream pressures, the input resistance and inertia, and the calculated 
outputs of gain and phase (Fig. C-3) . The upstream and downstream pressures, 
as well as the resistances and inertances, are tabulated in this manner for 
each flow to simplify checking for proper input of the data. For each pres- 
sure mode, the printout shows the flows in and out, the volume and acoustic 
velocity, and the calculated outputs of gain and phase (Fig. C-4) . Again, 
this format simplifies model checkout. 

Also shown in Fig. C-3 are the four model summations of the injector flow- 
rates in both lb/sec/psi and percent flow per percent chamber pressure. 


The output of magnitude and phase gives information on the amount of response 
throughout the system. For example, nodes 49 through 60 represent the pres- 
sures in the outer ring groove. For this case, the amplitude of pressure node 
51 is 0.9399 at a phase angle of 335.3 degrees which is nearly the same as 
the maximum chamber pressure input. Therefore, the volumes in the injector 
manifolds are not so large as to filter out the pressure oscillations in the 
ring grooves. This means that the injector flow oscillations can be affected 
by changes upstream of the injector face. If there had been little or no pres- 
sure response in the ring groove, the only feed system modifications that could 
alter the system would be changes to the injector orifices. Another comparison 
that lends insight into the system response is the amount of injector' flow to 
ring groove flow. The maximum amplitude for flow through the injector for the 

outer ring groove is 4.716 x 10 ^ lb/sec per psi while the maximum ring groove 
-4 

flow is 3.396 x 10 . Since the ring groove flow is significant compared to 

the injector flow, dams in the ring groove could significantly affect the 
ring groove pressure response, and therefore the injector flow oscillations. 
This type of analysis for a specific injector has to be performed for each of 
the ring grooves, and throughout the rest of the injector flow passages. 
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Before any conclusions can be reached concerning injector hydraulic coupling 
in a specific injector, both the oxidizer and fuel feed systems must be ana- 
lyzed. The relative injector flow oscillations between the oxidizer and fuel 
sides may indicate the likelihood of one side being the controlling side, or 
indicate if both sides are important. Injector modifications could then be 
analyzed for possible beneficial effects. 



SECTION IV 


MODEL CORRELATION 


Utilizing the model whose development is documented in Section III, three 
different test cases (injectors) for which test data exists and for which 
coupling between the injector hydraulics was known to occur were analyzed. 

The test cases selected were those cases for which solutions to the coupling 
problems were obtained and whose solutions could be analyzed using the 
developed computer model to post-predict the influence of the solution rela- 
tive to the potential for coupling. These test cases were: the Rocketdyne 

Space Shuttle OME Technology Injector, the Rocketdyne Lance XPL Booster 
injector, and the Aerojet Space Shuttle OME Technology injector. 

REVIEW OF TEST CASES 

Compilation of test histories for the three above mentioned injectors are 
summarized below. These injectors most clearly exhibit a "high-frequency" 
coupling between the injector hydraulics and the chamber dynamics. Injector 
manifolding details for these injectors have been previously discussed in 
Section II of this report and reference is frequently made to figures pre- 
sented earlier in Section II. 

Rocketdyne Space Shuttle OME Technology Injector 

The Rocketdyne SS/OME technology engine produces 6000 pounds thrust at a cham- 
ber pressure of 125 psia using N 2 0^/MMH propellants. Chamber diameter and con- 
traction ratio are 8.2 inches and 2.0, respectively. A like-doublet injection 
element pattern is utilized. The injector has a total of 186 primary elements 
(744 orifices) having orifice sizes of 0.032 to 0.038-inch diameter (ox) and 
0.028 to 0.033-inch diameter (fuel). Injector AP's are 56 psid (ox) and 62 
psid (fuel) . 

Details of the Rocketdyne OME technology injector were shown earlier in Figs. 
10-14. The injection orifices are fed by a total of 10 ring grooves (5 ox and 
5 fuel) behind the 8.2-inch-diameter injector face. The oxidizer ring grooves 
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are fed through slanted feeder passages from a central oxidizer manifold. The 
fuel ring grooves are fed through slanted feeder passages from an annular fuel 
manifold. 

A dual mode quarter-wave absorber with a contoured entrance (see Fig. 10, 

■> Acoustic Cavity) was included in the original OME technology injector design. 

The annular slot was partitioned radially by 12 dams. Eight of the twelve 
individual slots were "tuned" to damp the IT acoustic mode of the chamber 
t having a 1.60-inch effective slot depth and 0.5-inch slot width (14.8 

percent open area). The remaining 4 slots were tuned to damp the 3T, 1R 
acoustic modes of the chamber and had a 0.78-inch effective slot depth 
and 0,5-inch slot width (7.4 percent open area). With this absorber con- 
figuration, the OME combustor was stable until the combustion was perturbed 
by an explosive bomb. A 25-percent occurrence of dynamic instability was 
obtained, a 2600 Hz mode of instability being recorded. Analysis of this 
instability mode* indicated a coupling between the injector hydraulics and 
the combustion process (Ref, 13). The addition of three radial dams in 
the annular fuel manifold succeeded in achieving a dynamically stable 
combustor, 

A slight alteration in slot width to obtain open areas of 12 percent (IT) 
and 6 percent (3T, 1R) , however, resulted again in a 12-percent occurrence 
of dynamic instability. The 2800 Hz instability mode was again believed 
indicative of hydraulic coupling and the addition of dams in the injector 
ring grooves was considered. No further testing of this injector was done, 
however . 


*The mode appeared to have the angular pressure distribution of a first 
tangential mode with the nodal position influenced by location of the 
fuel inlet. 
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Rocketdyne Lance XRL Booster Injector 

The Lance XRL booster engine produces 42,000-lb-thrust at a chamber pressure 
of 950 psia using IRFNA/UDMH propellants. The annular combustor has a 12.5 
inch O.D. and a 7.0-inch I.D. Contraction ratio is 2.8. The injector pattern 
consists of four concentric rows of unlike doublets that impinge 0.250-inch 
from the injector face. The orifice hole sizes are 0.073-inch diameter for 
the oxidizer and 0.0515 inch for the fuel. Injector AP's are 184 psid (ox) 
-and 140 psid (fuel) . Figure 7 shows the frontside view of this injector 
face. Top and side views of the overall injector are provided by Figs. 8 
and 9. 

During development of the flightweight engine system, oscillations with fre- 
quencies from 1000 to 9000 Hz were incurred. Theoretical and empirical 
studies were conducted which culminated in the successful stabilization of 
this engine (Refs. 4 and 14). 

A four-vane, 5 . 75- inch-long injector face baffle (chordal alignment) was in- 
cluded in the injector design to preclude occurrences of tangential insta- 
bility. A combination of "L" shaped Helmholtz resonators was included in the 
injector design as shown in Fig. 26 to preclude oscillatory frequencies be- 
lieved associated with second-baffle compartment modes (4400 Hz), the fourth 
tangential model (4800 Hz), the sixth and eighth tangential mode, and the 
first radial mode, and thus provide dynamic stability of the flightweight 
engine. The entrance to the modified Helmholtz cavity is located at the in- 
jector face and the cavity volume is located in the ablative thrust chamber 
body. The XRL booster engine incorporates 16 modified Helmholtz resonators 

(four per baffle compartment) tuned to 6120 Hz. Each resonator has an aper- 

2 

ture area of 0.238 in. , an aperture length of 0.100 inch, and a cavity vol- 

3 

ume of 0.552 in. . The absorber open area is thus 6.5 percent of the total 
injector area. The acoustic absorbers are designed to provide strong atten- 
uation over a wide band on either side of the tuned frequency. The absorbers 
proved completely successful in damping the modes of instability for which 
they were designed. 

Despite the presence of both baffles and acoustic absorbers, the engine re- 
peatedly exhibited a unique 1300 Hz, sinusoidal, linear instability. The 
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Figure 26. Injector Fnd Absorber 

observed instability was determined from test data and model simulation to 
be a chamber resonant mode with flow under the baffles combined with hydraulic 
coupling of the injector (Fig. 27 ). For this mode, pressures in adjacent 
baffle compartments were 180 degrees out of phase (observed from pressure 
measurements in the various compartments) . This resonant chamber mode is 
a slight variation of the second tangential mode in that, at the injector 
face, there is no flow but near the chamber throat the flow is tangential. 

This mode has a frequency lower than the second tangential (2490 Hz) due to 
the additional flow path length to get under the baffles. 

The oxidizer and fuel injector feed systems had four individual feeder arms 
that distributed propellants to ring grooves. From the ring groove, the pro- 
pellant flowed through the injector orifices into the combustion chamber. A 
"lumped parameter" injector model was made for each of the propellant mani- 
folds in the Lance XKL booster (Fig, 28) to analyze the unique 1300 Hz mode 
of instability. The model employed a numerical method of solution of the 
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Figure 27. Schematic of XRL Booster Injector Without Dams Showing 
Ring Groove Crossflow Model Causing Coupling With 
1300 Hz Anomalous Second Tangential Mode 
















continuity and momentum equations to obtain a predicted injector response to 
a pressure oscillation mode imposed at the injector face. 

Thus, the response of the various elements was affected by tangential 
flow in the ring grooves. The ring groove pressures tended to follow 
the oscillations in chamber pressure. Therefore, in the ring grooves 
just above a baffle, the pressures tended to be out of phase which 
produced tangential ring groove flow. The net effect was to produce 
the pressure distribution shown in Fig, 29 . The shaded area then re- 
presents the oscillatory injector pressure drop which produced injector 
flow oscillations. By putting in ring groove dams just above the location 
of the baffles (location of dams above the baffles is very important), the 
pressure distribution shown in Fig. 30 was predicted by the model. 

The shaded area representing the oscillation in injector pressure drop 
was drastically reduced which made injector flow nearly constant. In- 
corporation of the ring groove dams in the XRL engine completely eliminated 
the 1300 Hz mode of instability with no change in injector pressure drops 
or other engine operating parameters. The Lance XRL booster engine is in 
production today and over 1000 have been produced, operating stably with 
ring groove dams, acoustic absorbers, and baffles. 

Aerojet Space Shuttle OME Technology Injector 

The Aerojet engine produces 6000 Ibf thrust at a chanter pressure of 
125 psia, using nitrogen tetroxide, N^O^, and monome thy 1/hydrazine, MMH 
propellants. Details of the Aerojet injector are shown in Figs. 15 
and 16. An 867-element X-doublet platelet pattern (see Fig. 31) on the 
8.2-inch-diameter injector face is fed by a total of 15 ring grooves (8 
fuel and 7 ox) . An outer ring circumferential fuel manifold feeds an 
inner ring fuel manifold through 24 drilled holes. Fuel flows from the 
inner fuel ring through 3 downcomer slots to 3 pie-shaped manifolds. Down- 
comer slots from the pie manifolds feed each fuel ring. The oxidizer is 
fed through a central manifold to 3 pie-shaped manifolds. As in the case 
of the fuel, downcomer slots from the ox pie manifold feed each ox ring. 

A baseline injector-acoustic cavity configuration was established, which 
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Figure 30, Response With Ring Groove Dams 






consisted of an X-doublet platelet injector, in conjunction with a circumfer- 
ential cavity housing. The cavity housing holds eight 1-T (nominal 18% open 
area) and four 3-T acoustic cavities (nominal 9% open area) . Cavities were of 
the quarterwave slot type. 

The primary mode of instability with the X-doublet pattern was resurging, which 
consists of periodic burst of a spinning 1-T mode at about 2600 Hz. The fre- 
quency of these bursts is about 400 Hz. Three changes to the OME engine were 
tested to determine their effect on eliminating resurging. Two of these changes 
dealt with the OME injector, while the other dealt with the combustion chamber 
acoustic cavities. 

The most successful of the three changes was found to be the installation of 
ring groove dams. These dams were installed at the three null points in each of 
the outer 12 ring grooves. The insertion of a ring in the fuel circuit torus as 
shown in Fig. 32 to increase the injector pressure drop was another change incor- 
porated in the injector to eliminate resurging. However, this had no demonstrated 
effect in changing resurge behavior. Finally, acoustic cavity geometry was al- 
tered considerably as the last of these changes. Both cavity inlet area and over- 
lap (see Fig. 33) were altered. The results showed that decreasing inlet area 
tended to give more acoustic mode instabilities with less resurging, but that 
increasing overlap improves stability. For more detail on the Aerojet OME in- 
jector, see Ref. 8. 

Thus, like the Rocketdyne OME technology injector, the Aerojet injector also 
exhibited a 2600 Hz instability whose mode was feed system coupled. Instabil- 
ity was highly resistant to suppression with acoustic cavities alone, but was 
influenced by the cavity entrance configuration. As in the Rocketdyne case, 
instability suppression was most successfully achieved through the use of dams , 

PREPARATION OF MODEL INPUT DATA 

Model input data was prepared for the three chosen test cases. This 
included computing the capacitances, inertances, and linearized flow 
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Figure 32. Fuel Circuit Pressure Drop Ring 
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Figure 33. Cavity Overlap Results 












resistances for the "lumped" fuel and oxidizer sides of each injector test 
case as described in the previous section. Schematics of these lumped in- 
jectors are shown in the discussions to follow. 

Rocketdyne Space Shuttle OME Technology Injector 

Modeling of this injector was concentrated on the analysis of the fuel side 
due to the fact that all injector fixes were associated with this side. 

The fuel side of this injector was "lumped out" as schematically shown in 
Fig. 34 , with 78 pressure nodes (shown by rectangles), 168 flow segments 
(shown by lines) and 39 pressure input locations (shown by circles). The 
black triangles in this figure and all further injector schematics denote 
the locations of injector dams which were later incorporated into the in- 
jector design after unstable operation occurred. The angular locations of 
these pressure nodes, flow segments, and dams can be determined from Fig. 35. 
By assuming negligible pressure drops in flow segments 100-138 for the non- 
dammed case, the fuel side was reduced to a simple network of series and 
parallel flows as shown in Fig. 36 . This assumption is made from the fact 
that the time-averaged steady-state flow is very small in these segments 
giving a time-averaged pressure drop of zero. Based on a total flowrate of 
7.27 lbm/sec an overall injector pressure drop of 51 psid was calculated. 

This compares to a reported pressure drop of 62 psid which gives a 17% error. 

The oxidizer side of the Rocketdyne OME technology injector was described 
as shown in Figs. 37 and 38 . The oxidizer side consists of 69 pressure 
nodes, 129 flow segments, and 31 pressure input locations. By assuming 
negligible pressure drops in flow segments 70-99 the oxidizer side was 
reduced to a simple network of series and parallel flows as illusrated 
in Fig. 39 . Based on a total oxidizer flowrate of 12 lbm/sec, an 
overall .injector pressure drop of 47 psid was calculated. This compares to 
a reported pressure drop of 56 psid for a 17% error. 

Rocketdyne Lance XRL Booster Injector 

The fuel side, schematically shown in Fig. 40, consists of 53 pressure 
nodes, 88 flow segments, and 28 pressure input locations. The angular 
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Figure 34 


Lumped Parameter Representation of Fuel Side of Roclcetdyne OME Technology 
Injector Showing Flows Between Pressure Nodes 



































































































Annular Fuel Manifold (9 lumps) 


.Ring Grooves No. 8, 10 (9 lumps apiece) 


Ring Grooves No. 4, 6 (6 lumps apiece) 



Location of 
Manifold Dams 
(typical, 120° 
spacing) 


Ring Groove No. 2 (3 lumps) 


Figure 35. Distribution of Lumps for Fuel Side of 
Rocketdyne OME Technology Injector 
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Figure 36. Simplified Resistance Network for Fuel Side 


Fuel Reservoir 



Rocketdyne OME Technology Injector 







































































Figure 38. Distribution of Lumps for Ox Side of 
* Rocketdyne OME Technology Injector 
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Figure 40. Lumped Parameter Representation of Fuel Side of Rocketdyne 
XRL Injector. Showing Flows Between Pressure Nodes 




















locations of these lumped constituents can be determined from Fig. 41. 

By assuming negligible pressure drops in flow segments 17-44 and 77-88 the 
fuel side was reduced to the simple network as shown in Fig. 42 . Based on 
a total fuel flowrate of 33.3 lbm/sec an overall injector pressure drop of 
129 psid was calculated. This compares quite well to a reported pressure 
drop of 140 psid which yields an error of 8%. 

The oxidizer side, schematically shown in Fig. 43 , consists of 61 pressure 
nodes, 108 flow segments, and 36 pressure input locations. The angular 
placement of these constituents can be determined from Fig. 44 . By assuming 
negligible pressure drops in flow segments 37-72 the oxidizer side was reduced 
to a simple network of series and parallel flows as illustrated in Fig. 45. 
Based on a total oxidizer flowrate of 133 lbm/sec, an overall injector pres- 
sure drop of 211psid was calculated. This compares to a reported pressure 
drop of 184 psid which gives a 15% error. 

Aerojet Space Shuttle OME Technology Injector 

The fuel side, schematically shown in Fig. 20 , consists of 99 pressure nodes, 
223 flowrates, and 52 pressure input locations. The angular locations of 
these lumped constituents can be determined from Fig. 19. By assuming negli- 
gible pressure drops in flow segments 2-16, 32-46, 53-73, and 122-172, the 
fuel side was reduced to the simplified network shown in Fig. 24. An overall 
fuel side pressure drop of 58 psid was calculated based on a total fuel flow- 
rate of 7.19 lbm/sec. This compared very well with a 57 psid pressure drop 
reported by Aerojet which gives an error of 2% between the two values. 

The oxidizer side of the Aerojet OME technology injector was schematically 
described as shown in Fig. 46. The oxidizer side consists of 71 pressure 
nodes, 171 flow segments, and 50 pressure input locations. The angular 
position of these components can be determined from Fig. 47 . By assuming 
negligible pressure drops in flow segments 11-31 and 75-122, the oxidizer 
side was reduced to a simple network of series and parallel flows as illus- 
trated in Fig. 48. Based on a total oxidizer flowrate of 11.86 lbm/sec an 
overall injector pressure drop of 56 psid was calculated. Thxs compares 
to a reported pressure drop by Aerojet of 44 psid for a 27% error. 
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Figure 41. Distribution of Lumps for Fuel Side of 
Rocketdyne Lance XRL Booster Injector 
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Figure 43. Lumped Parameter Representation of Oxidizer Side of Rocketdyne 
XRL Injector Showing Flows Between Pressure Nodes 



/I 


u™ 

-io * 

f 

11 


*■ 

aa 

1 [ 

J5 | 


*7 

>1 

3 J.2L3Z 
- . - A 

£j 33 -5^ 5V 


( 

) ( 

i <£ 

& & a 

■ ■ SJ 
1 









V 

1 

I 

. 

r 




/#5 /•¥ 


ITT] 

Lad 


l!LJ 

rn 

, s i 


(35) jy) 


Manifolds 


3 7 _e. 3 g *u 3 ? -S. yo -!z_ 3 Ring Groove 1 

Oil 

j9*. 

T s , LeU „ L_l> DRing Groove 3 

© © 



nJTsi LziJ LO ±£ D Ring Groove 5 

A St 1 lz tz 

fas) fa?) 






















































Figure 45, Simplified Oxidizer Side 



XRL Injector 



Figure 46. Lumped Parameter Representation of Oxidizer Side of Aerojet OME Technology 
Injector Showing Flows Between Pressure Nodes 
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Figure 47. Distribution of Lumps for Oxidizer Side 
of Aerojet OME Technology Injector 
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Figure 48. Simplified Oxidizer Side, Aerojet 
OME Technology Injector 




CORRELATION OF MODEL 


The engine hydraulic stability model was run for the Rocketdyne OME technology 
and XRL engines and for the Aerojet OME technology engine. Complete sets of 
input and output data are shown for each side of the three correlation injec- 
tors in Appendices C through H.* 

Rocketdyne Space Shuttle OME Technology Injector 

The Rocketdyne OME technology injector exhibited a IT mode of instability on 
25% of the tests where bombs were detonated in the chamber. By adding three 
radial dams in the annular fuel manifold, the instability was eliminated. 

The same injector with the dams and a different acoustic cavity design had a 
12% incidence of instability when bombed. This means that the incorporation 
of dams did reduce the loop gain, but the reduction was probably not a big 
reduction or the modification of the acoustic cavities would not have caused 
the instability to reappear. 

The lumped parameter representations for the fuel and oxidizer sides of the 
Rocketdyne OME technology injector are shown in Figs. 34 and 37. Input for 
the fuel side without the dams is shown in Appendix E. The resistance values 
for flows 130, 133, and 136 (see lumped parameter description shown in Fig. 34) 
were set equal to 10^ when simulating the addition of manifold dams to the 
fuel side of the injector. A standing first tangential mode (v = 2600 Hz) 
in the combustor was input to the injector face. A summary of the injector 
flow results is shown in Table III. This table has the four different total 
injector flow summations with the output in both lb/sec/ptii and % flow/% P c> 

The results presented in Table III show about a 17% reduction in fuel flow 
with the radial dams. Plots for the outer fuel ring groove with and without 
dams are shown in Figs. 49 and 50. Therefore, it is concluded that adding 
the annular dams reduces the hydraulic coupling with the combustion chamber 
thus increasing, to a limited extent, the stability of the combustor. 


* These appendices show results for the injectors of interest without the 
inclusion of dams or any other injector "fix." 
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TABLE III 



Total 

Injector 

Flow 



Lb /Sec 
Psi 

% Flow 
% P 

c 


Fuel Side 
No Dams 

.00633 

.109 


Fuel Side 
With Dams 

.00521 

.0896 


Ox Side 

; 




No Dams 


00551 


0574 
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Rocketdyne Lance XRL Booster Injector 

The lumped parameter representations for the fuel and oxidizer sides of 
the XRL are shown in Figs. 40 and 43. The XRL had a 1300 Hz mode with 
chamber pressure amplitudes of 200 to 300 psi peak-to-peak on every test. 
Incorporation of ring groove dams at four locations in each ring groove 
(directly above chamber baffle locations) eliminated the 1300 Hz mode. 

A summary of the results for the oxidizer and fuel sides are shown in 
Table IV. The oxidizer side amplitude was reduced to less than one-third 
of the original value while the fuel side w r as reduced to about 70 percent 
of its original value. Typical plots for the oxidizer and fuel side with 
and without dams are shown in Figs. 51 through 54. The model shows a sub- 
stantial reduction in hydraulic coupling which is consistent with the ob- 
served results. 

Aerojet OME Technology Injector 

The lumped parameter representations for the fuel and oxidizer sides of the 
Aerojet injector are shown in Figs. 20 and 46. The Aerojet injector had a 
spinning IT mode at about 2600 Hz with resurging at about 400 Hz. Installa- 
tion of dams at the three null points in each of the outer 12 ring grooves 
was most successful in eliminating the instability . 

Results of the model for the fuel and oxidizer side with and without dams 
is shown in Table V. These results show that the inclusion of dams on the 
oxidizer side had very little effect on the injector flowrates. However, 
inclusion of dams on the fuel side increased the injector oscillatory flow 
by a factor of 1.4. 

Initially, this seemed inconsistent with the observed results. However, a 
re-examination of the resurging type of instability experienced by Aerojet 
was made to determine whether the model prediction of increased injector 
coupling with a spinning IT chamber mode did indeed correlate with the ex- 
perimental results. The resurging instability is reported to manifest itself 
as periodic low frequency (~400 Hz) burst of high frequency (>2000 Hz) 
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TABLE IV. XRL INJECTOR RESULTS 



XRL Fuel Side 
No Dams 


XRL Fuel Side 
With Dams 


XRL Ox Side 
No Dams 


XRL Ox Side 
With Dams 


Total 
Inj ector 
Flow 

Total Vector 
Injector 
Flow 

Total Injector 
Flow Proportioned 
by PC Amplitudes 

Total Vector 
Injector Flow 
Proportioned by 
PC Amplitudes 

Lb /Sec 

% Flow 

Lb /Sec 

% Flow 

Lb/ Sec 

% Flow 

Lb /Sec 

% Flow 

Psi 

% P 

Psi 

% P 

Psi 

% P 

Psi 

% P 


c 


c 




c 

.0551 

1.34 

.0544 

1.32 

.0551 

1.34 

.0544 

1.32 

.0397 

.966 

.0394 

.958 

.0397 

.966 

.0394 

.958 

.1246 

.911 

.1242 

.907 

.1246 

.911 

.1242 

.907 

.0386 

.282 

.0270 

.198 

.0386 

.282 

.0270 

.198 
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XRL Injector - No Dams 
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Figure 52. Results for Fuel Ring of 
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Total 
Inj ector 
FZ :<t 



Lb/ Sec 
Psi 

% Flow 
% P 

c 

Fuel Side 
No Dams 
Spinning IT 

.0194 

,/ .337 

Fuel Side 
With Dams 
Spinning IT 

.0267 

.464 

Ox Side 
No Dams 
Spinning IT 

.0636 

.670 

Ox Side 
With Dams 
Spinning IT 

.0629 

.663 


,T OME TECHNOLOGY INJECTOR RESULTS 


Total Vector 
Inj ector 
Flow 


Lb /Sec 


Psi 


.0167 


Total Injector 
Flow Proportioned 
by PC Amplitudes 


Total Vector 
Injector Flow 
Proportioned by 
PC Amplitudes 


% Flow 

Lb /Sec 

% Flow 

Lb /Sec 

% P 

c 

Psi 

% P 

c 

Psi 

.290 

.0149 

.259 

.0130 

.396 

.0211 

.367 

.0179 

.619 

.0545 

.574 

.0504 

.616 

.0532 

.561 

.0496 




















































instability. The resurge is reported in Ref. 9 to originate with a low 
amplitude spinning IT wave which grows rapidly in amplitude and frequency 
until it becomes a detonation wave. The detonation wave is made possible 
because of an accumulation of unburned propellant (due to poor mixing in- 
herent with the x-doublet injector pattern) within the chamber near the 
injector face. The detonation wave makes one circuit of the chamber wiping 
out all propellant in its path. Pressure then decays (with acoustic ringing) 
and steady-state flow is re-established. Acoustic cavities were found to be 
unable to damp out the resurge and, in fact, were believed to foster the 
resurge by virtue of damping the high frequency modes which would normally 
promote mixing and combustion of the unburned propellant. Indeed, it is 
reported that when IT damping -was removed from the system*, the resurging 
disappeared to be replaced by a normal (standing) IT or 2T acoustic mode. 
According to the engine hydraulic stability model, the addition of ring 
groove dams increases (not decreases) injector coupling with the spinning 
IT chamber mode. This is interpreted as decreasing the IT damping and 
promoting the mixing and combustion of the unburned propellant which other- 
wise would accumulate near the injector face and ultimately contribute to 
the resurging instability.** Therefore, it is concluded that the model 
does correlate with the observed test data. 

Typical plots for the outer ring grooves with and without dams are shown in 
Figs. 55 through 58, 

As described above, the engine hydraulic stability computer model was run for 
each of the "correlation" injectors with both pre-fix and post-fix injector 
model input. Examination and analysis of the model output reveals that the 
computer model successfully predicted that the fixes applied to each correla- 
tion injector would increase combustor stability with respect to the instabi- 
lity mode actually observed. 

*By virtue of changes in the acoustic cavity. 

**Sustainment of the IT or 2T mode is precluded by the presence of the 
acoustic cavities . 
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Figure 56. Results for Outer Fuel Ring Groove of 
Aerojet Injector - Dams 
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SECTION V 


SENSITIVITY ANALYSIS 

A systematic parametric variation of some major injector design variables was 
made to gauge the effect of the individual parameters on injector gain and 
to identify the most sensitive variables. The objective of this task was to 
develop, to the maximum extent possible, generalized design criteria to pre- 
clude coupling between the injector hydraulics and the combustion process. 

The "sensitive variables" whose effect was investigated (using the computer 
model) included: 

1 . Propellant temperature 

2 . Frequench of a particular mode 

3. Resistance values of the injector orifices 

4. Inertance values of the injector orifices 

5. Ring groove area (i.e., capacitance and inertance of the ring 
groove nodes) 

6. Chamber instability mode 

7. Presence of dams 

The effect of these variables was computed for the fuel and oxidizer sides 
of both the Aerojet and Rocketdyne OME technology injectors. Model predic- 
tions are summarized in Tables VI through IX. These same predictions are 
also displayed graphically in Figs. 59 through 67 where the gain plotted is 
based on the total vector injector flow proportioned by P^ amplitude. 

It is evident from Tables VI through IX and Figs. 60, 62, 64, and 66 that 
any error in linearized resistance input to the model due to disagreement 
between calculated and experimental injector AP* is unimportant from a 


* See Section IV. 
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TABLE VI 


MODEL PREDICTION OF AEROJET FUEL SIDE GAIN* 


H 

H* 

to 




o 

nS 

SC £3 

‘ & 


Off Nominal Condition** 

Ring Groove Dana 


Ho 

No 

No 

Yea 

Yea 

Yea 

Yes 

Variable 

Input Value 

No 



Tot. InJ. Flow 

Tot. Vet InJ. Flow 

Tot, a P £ Anp. 


Tot InJ. Flea 

Tot Vet InJ. Flow 

Tot o P £ Anp. 

Tot Vet a P Aap. 
c 

Noal 

nal** 

.0194 

.0167 

.0149 

,0130 

.0267 

,0228 

.0211 

.0179 

InJ. temp, 

70°F 

.0129 

,00637 

,0104 

,00637 

,0160 

.00414 

,0127 

,00424 

InJ. temp. 

200°F 

.0151 

.0113 

,0112 

, 00863 

,0231 

,0183 

.0181 

.0143 

X»J. temp. 

300°F 

.0265 

,0242 

,0210 

,0192 

,0308 

.0278 

.0246 

.0221 

Frequency 

2200 Hr 

.0117 

.00304 

,00866 

.00315 

,0190 

.0113 

.0149 

.00895 

Frequency 

3000 Hz 

.0263 

.0261 

.0224 

,0207 

.0325 

.0293 

.0259 

.0233 

Resistance 

0.75 R noa 

.0207 

.0167 

,015$ 

,0127 

,0296 

.0238 

.0232 

.0185 

Resistance 

1.25 R noa 

,0184 

,0164 

.0143 

,0129 

.0243 

.0214 

.0193 

.0169 

Orifice lnertance 

0.75 I noa 

,0188 

.0161 

.0144 

,0126 

.0256 

.0218 

.0202 

.0171 

Orifice lnertance 

1.25 I noe 

.0200 

.0173 

,0154 

.0135 

.0280 

.0239 

.0221 

,0188 

Orifice lnertance 

2.0 1 noa 

.0224 

.0197 

,0173 

.0153 

.0324 

.0281 

.0257 

.0222 

Orifice lnertance 

3.0 I noa 

.0267 

.0241 

,0208 

,0189 

.0392 

.0348 

.0312 

.0276 

Ring groove area 

-25X noa*«* 

.0196 

.0169 

.0150 

.0131 

.0249 

.0214 

.0195 

.0167 

Ring groove area 

+25X noa**** 

.0190 

.0164 

.0147 

.0128 

.0284 

.0241 

.0227 

.0191 

Hod* 

Standing IT 

.0145 

.00938 

,00736 

.00654 

.0196 

• .0124 

.0101 

.00900 

Mode 

Standing 2T 

.0185 

.0148 

,00958 

.00810 

.0181 

,0142 

.00964 

.00785 


the four single value ■ (natation a. 
z - 2600 He, spinning IT mode 


•Gain, (lb/aec) /pal , la ahown for aaah of 

••Nominal conditions Include T. . » 230pr, 
lnj 

•••0.75 C 1.333 I 

non, not 
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TABLE VII. MODEL PREDICTION OF AEROJET OX SIDE GAIN* 


Off IfoBlnal 

Condition** 

Ring Croove Dast 

Varlahle 

Input Value 

No 

No 

No 

No 

.. ......... 

Yes 

Yea 

Yea 

. 

Yea 



Tot. InJ. Flow 

Tot. Vec. InJ. Flow 

Tot. a P c A np. 

Tot. Vec. a fmp. 

Tot. InJ. Flow 

Tot. Vec. InJ. Flow 

Tot. a Aarp 

Tot. Vec. a P c A*p. 

Noulna 


.0636 

.0587 

.0545 

.050, 

.0629 

.0584 

.0532 

.0496 

Frequency 

2200 Hr 

,0560 

.0512 

.0495 

.0441 

.0618 

,0542 

.0521 

.0462 

Frequency 

3000 He 

.0560 

.0522 

.0479 

.0446 

.0543 

.0516 

.0459 

.0438 

Resistance 

0.75 R no* 

.0829 

*0762 

.0710 - 

.0656 

.0827 

.0759 

.0701 

.0647 

Resistance 

1.25 R no* 

.0521 

.0481 

.0446 

.0411 

.0514 

.0481 

.0434 

.0407 

Orifice lnertance 

0.75 I no* 

.0627 

.0574 

.0536 

,0492 

.0622 

.0574 

,0526 

.0487 

Orifice lnertance 

1.25 I non 

.0629 

.0585 

,0539 

.0501 

.0618 

.0577 

.0524 

.0491 

Orifice lnertance 

2.0 X no* 

.0533 

,0504 

.0457 

.0432 

.0518 

.0492 

.0440 

,0418 

Orifice lnertance 

3.0 I no* 

.0390 

.0374 

.0334 

.0320 

.0378 

.0364 

,0321 

.0310 

Ring groove area 

-25X no* 

,0609 

.0563 

.0519 

.0481 

.0598 

.0560 

.0509 

.0477 

Ring groove area 

+25X no* 

.0661 

.0610 

,0568 

.0524 

.0653 

.0603 

.0559 

.0517 

Mode 

Standing IT 

.0432 

.0388 

.0269 

.0256 

.0429 

.0389 

.0264 

.0253 

Mod* 

Standing 2T 

,0219 

.0164 

.0127 

.00970 

.0218 

.0134 

.0140 

.00844 


•Gain (lb/aec)/pai, la shown for each of the four tingle -value a taxations 
**No*inal conditions Include f * 2600 He, spinning IT *odet 
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TABLE VIII. MODEL PREDICTION OF ROCKETDYNE FUEL SIDE GAIN* 


Off Nominal Condition** 

Annular Manifold Dams 

Variable 

Input Value 

No 

No 

No 

No 

Yea 

Yes 

Yes 

Yes 



Tot. Inj, Flow 

Tot. Vec. Inj Flow 

Tot, a P c Asp. 

Tot. Vec. a P c Amp. 

Tot. Inj. Flow 

Tot. Vec. Inj. Flow 

Tot- a P c Amp. 

Tot. Vec. a P c Amp, 

Nominal 

. 

.00633 

.00632 

,00442 

.00441 

.00521 

.00508 

.00371 

.00367 

Frequency 

2200 Hr 

.00818 

.00817 

.00567 

.00566 

,00734 

.00731 

.00517 

.00516 

Frequency 

3000 Hr 

.00460 

,00457 

.00325 

.00324 

.09331 

.00224 

.00212 

.00185 

Resistance 

0.75 R nom 

.00637 

.00636 

.00445 

.00445 

.00523 

.00511 

.00373 

.00370 

Resistance 

1.25 R non 

,00627 

.00626 

.00437 

.00437 

,00518 

.00504 

.00368 

.00364 

Orifice inertance 

0.333 I non 

.00856 

,00849 

.00601 

.00598 

.00715 

.00648 

.00497 

.00480 

Orifice inertance 

0.50 I nom 

.00786 

.00782 

.00551 

.00550 

.00652 

.00606 

.00457 

.00445 

Orifice inertance 

.0,75 I non 

. .00701 

.00699 

.00490 

.00490 

.00578 - 

.00552 

.00409 

.00402 

Orifice Inertance 

1.25 I non 

.00577 

.00576 

,00402 

.00401 

.00477 

.00470 

.00340 

.00338 

Ring groove area 

-251 non 

.00623 

.00622 

.00437. 

'.00437 

.00512 

.00494 

.00366 

.00361 

Ring groove area 

4-25T non 

,00661 

,00659 

.00463 

.00463 

.00520 

.00375 

.00375 

.00373 

Mode 

Spinning IT 

.00988 

.00986 

, 00863 

.00882 

.00827 

.00805 

.00750 

.00735 

Ring groove dams 

- 

,00524 

.00510 

.00382 

.00378 

,00468 

.00318 

.00310 

.00266 


•Gain (lb/aec)/pai, is shcum for each of the four single value emanations, 

••Nonlnal conditions include T, . ■ 200°F, f • 2600 Hr. standing IT node, 
inj 
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TABLE IX. MODEL PREDICTION OF ROCKETDYNE OX SIDE GAIN* 




•Gain (lb/aec)/pai, is shown Cor each of Che four single value simatlone. 
••Nominal conditions Include f ■ 2600 Hz, standing IT mode. 























Gain, (]b/sec)/psl 














Gain, (lb/sec)/ps! Gain, (Ib/se 



0 1 2 3 



Ring Groove Area 


Figure 61, Effect of Orifice Inertance and Ring Groove Area on 

Gain for Fuel Side of Aerojet OME Technology Injector 




Gain, (1b/sec)/psi Gain, (lb/sec)/psi 



2200 2600 3000 


no dams 

ring groove dams 



0.75 1-0 1 .25 


R/R 


nom 


I 


Figure 62 • Effect of Frequency and Resistance on Gain for 

Oxidizer Side of Aerojet OME Technology Injector 
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Ring Groove Area 


Figure 63 . Effect of Orifice Inertance and Ring Groove Area on 

Gain for Oxidizer Side of Aerojet OME Technology Injector 





















Ring Groove Area 


Figure 65. Effect of Orifice Inertance and Ring Groove Area on 

Gain for Fuel Side of Rocketdyne OME Technology Injector 



























model correlation standpoint. That is, the difference in gain between the 
base case and the base case with injector fix is not a function of the linear- 
ized resistance over the range of uncertainty in injector AP. 


It is observed in Figs. 59 through 67 that the general effect a specific para- 
meter has on the gain varies between the Aerojet and Rocketdyne injectors 
and, indeed, between the fuel and oxidizer side of the Aerojet injector itself. 
In order to gain an understanding of why the predicted gain varies in a par- 
ticular manner for various changes, an evaluation of the expected system 
responses was made. It is recalled that the linearized equation for flow- 
rate across a segment is: 


w 



R + IS 


where 


P u is the oscillatory upstream pressure 

P^ is the oscillatory downstream pressure 

R is the linearized resistance 

I is the inertance 


(33) 


Assuming that the flow segment is judt upstream of the combustion chamber 

(the orifice) P , = P , and setting S = jco, Equation (33) may also be written 
d c 

as 


w = 1/R = 1/R 

P -P 1 + (I/R) jco 1 + jWw b 


(34) 


where 

to is frequency in radians/sec 

co. is termed the break frequency = R/l 
b 
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Figure 68 shows a plot of log gain* versus log frequency. At low frequencies, 
i.e., to much less than to b> the gain &/ (P^-P ) becomes simply 1/R since 



Log Frequency (to) 


Figure 68. Gain versus Frequency 


5> Xn Fig. 68 gain is denoted by the expression w/0? ) whereas in the computer 

model output (see Tables VI through IX) gain is Henoted by £/P c . The latter 
term is employed in the model since it most effectively represents that quantity 
by which an injectors sensitivity to coupled instability is expressed. The 
term P u in the gain shown in Fig. 68 is, of course, an additional variable 
whose value (for fixed ? ) is dependent upon R, I, and to. 


the ico/co, term in the denominator of Eq. (34) becomes small compared to 
b 

unity. At high frequencies, i.e., w much greater than the gain becomes 
a function of both the inertance and the frequency (1/Iw) since be- 

comes large compared to unity. For frequencies near the break frequency, 
resistance, inertance, and frequency all affect the gain. 

-The break frequencies for the fuel and oxidizer sides of both the Aerojet 
and Rocketdyne OME technology injectors (calculated from injector orifice 
resistances and inertances) are: 


OME Technology Injector 

Break Frequency J 

Radians/ sec 

Hz 

Aerojet fuel side 

38,696 

6,159 

Aerojet Ox side 

23,500 

3,740 

Rocketdyne fuel side 

7,581 

1,206 

| Rocketdyne Ox side 

4,727 

752 


Looking at the sensitivity analysis results for the Rocketdyne injector 
shown in Tables XIII and IX, it is observed that gain is independent of 
orifice resistance, but a function of both orifice inertance and frequency. 
Since the results shown in Tables XIII and IX were obtained for a IT mode 
frequency of 2600 Hz, these observed model results are consistent with the 
anticipated results based on the consideration of the break frequency. That 
is, since the break frequencies for the Rocketdyne injector (1206 Hz- 
fuel and 752 Hz-ox) are significantly less than the model input frequency 
(2600 Hz), the flowrate gain should be primarily affected by inertance and 
frequency (higher inertance or frequency should reduce the gain) and not 
affected by resistance changes. Figures 64 - 67 show these trends. 

For the Aerojet injector, the break frequencies (6159 Hz-fuel and 3740 Hz-ox) 
are higher than the input frequency (2600 Hz) . Therefore, for a fixed input 
(P -P ) amplitude at 2600 Hz, it is expected that resistance increases will 
decrease the flow gain while inertance and frequency variations should have 
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little or no effect. The actual model results shown in Tables VI and VII do 
not agree well with these expectations. Thus, a closer examination of the 
model results was made. The model results showed that the injector pressure 
(]? ) was changing substantially, therefore overcoming the effect on flow gain 
due to resistance variation alone.* The model results also show (see Figs. 
60-63) that orifice inertance (I) and frequency (to) affect the flow gain as 
well. These observed results indicate that some other mechanism is affecting 
the model output for the Aerojet injector. 

Prior to proceeding with the identification of the mechanism believed res- 
ponsible for the unexpected trends in the gain of the Aerojet injector, an 
expedient effort was undertaken to. clarify the difference in gain between 
the Aerojet and Rocketdyne OME technology injectors. A total of seven addi- 
tional computer runs were made. These are summarized in Table X. Base runs 
for the above two injectors were made using identical input modes (spinning 
IT) and frequency (2600 Hz). Also, a common fuel injection temperature 
(230°F) was employed. Results from these four runs are shown in the first 
four lines of Table X. 

Table X shows that fuel side gains for the Aerojet injector are approximately 
two times as high as gains for the Rocketdyne injector. Oxidizer side gains 
for the Aerojet injector are approximately seven times as high as gains for 
the Rocketdyne injector. 

One possible explanation for these differences in gain is the difference 
between the ring groove volumes and orifice iner tances/resistances for the 
two injectors as summarized below in Table XI. 


*From Eq. (34), for to « to, , w / (P-P ) = 1/R. Therefore, 

^ a/ D u. C 

w/p c = K p u > Vr) . 
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TABLE X. COMPARISON BETWEEN MODEL OUTPUT FOR AEROJET 
AND ROCKETDYNE OME TECHNOLOGY INJECTORS 


1 njector 

Fuel or 
Ox Side 

Freq . 
Hz 

Mode 

■ 

Mod. 

Total 
l.nj. Flow 
(lb/sec)/ps! 

Total Vector 
InJ. Flow 
(lb/sec)/psl 

Tot. Inj. Flow 
Proport ioned 
by P c Ampl . 

( 1 b/sec)/ps i 

Tot. Vector 
Inj. Flow 
Proportioned 
by P c Ampl . 

(lb/sec)/psi 

Aerojet 

Fuel 

26 


i 



.0194 

.0167 

.01*49 

.0130 

Aerojet 

Ox 






.0636 

.0587 

.05*45 

.050*4 

Rocketdyne 

Fuel 






. 009^1 

.00938- 

.008*4*4 

.008*43 

Rocketdyne 

Ox 






.00880 

.00880 

.00739 

.00739 

Rocketdyne 

Fuel 





Aerojet 
Ring Grooves 
5 Or! flees 

.0102 

.00963 

.00912 

.00888 

Rocketdyne 

Ox 





.00970 

.00181 

.00813 

.00272 

Aerojet 

Ox 

1 

' 


f 

' 

Decreased 
Volume of 
P I e Man 1 fol d 
by Factor of 
10 

.OII 16 

.00305 

.0118 

.00287 

















TABLE XI. INJECTOR COMPARISONS 



Aeroj et 
Inj ector 

Rocketdyne 
Inj ector 

Fuel Ring Groove 
Volume, in. 3 

4.182 

1.80 

Oxidizer Ring^Groove 
Volume, in. 

5.39 

0.78 

Average Fuel Orifice 
Inertance 

2 . 3xl0~ 2 

6. 2x10" 2 

Average Oxidizer Orifice 
Inertance 

2.0xl0~ 2 

5 . 5xl0~ 2 

Average Fuel Orifice 
Resistance 

8.9xl0 2 

4.7xl0 2 

Average Oxidizer Orifice 
Resistance 

4 . 7xl0 2 

2.6xl0 2 


Two additional runs were therefore made for the Rocketdyne injector in which 
inputs representing the ring grooves and orifices were modified to simulate 
the substitution of the Aerojet injector ring grooves and orifices in the 
Rocketdyne injector. In general, this resulted in fuel and oxidizer side 
gains which fell in between the gains reported for the Aerojet and Rocketdyne 
injectors (see Table X). Although, this input modification shifted the results 
in the right direction (i.e., towards the Aerojet injector) the extent of the 
shift i^as not nearly as great as was expected (especially in the case of the 
oxidizer side). Consequently, the difference in predicted gain between the 
Aerojet and Rocketdyne injectors could not be explained totally from the stand- 
point of ring groove volume and orifice inertance/resistance differences. 

It was thus postulated that whatever mechanism was responsible for the un- 
expected trends of the Aerojet injector was also responsible for the qualita- 
tive difference in gain between the Aerojet and Rocketdyne injectors. 
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In order to determine what this mechanism was, a detailed evaluation of the 
injector inertances and volumes upstream of the orifices was made* Schematic- 
ally, the oxidizer and fuel sides of the Aerojet injector are as shown in Fig. 69 
The pie manifold feeders have a high enough inertance that they have little 
oscillatory flow in the frequency range of interest. The ring groove volumes 


Feeders (inertance) 


Pie Manifold (volume) 


Figure 69. 



Aerojet 


Downcomers (inertance) 

Ring Groove (volume) 

Injector Orifices (resistance 

inertance) 

Injector System Schematic 




are small compared to the pie manifold volumes. Therefore, this system can be 
closely approximated by a simplification (Fig. 70) that is essentially a 


Pie Manifold + Ring Groove, 
(volume) 


Downcomers + Injector (inertance) 


Figure 70. Aerojet Injector Simplified System Schematic 
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Helmholtz resonator. The resonant frequency of a Helmholtz resonator is: 


f 


- 5 _ I 

2 tt V g I 


where f 
a 
V 

g 

I 


resonant frequency, Hz 

acoustic velocity - in. /sec 

volume, in. ^ 

386.4 lb in . / lb c sec 2 
m f 2 

inertance, lb^ sec /lb in. 
5 f m 


Therefore, the calculated Helmholtz resonant frequencies for the Aerojet 
oxidizer and fuel sides are 2460 Hz and 4290 Hz, respectively. This means 
that peaks in the flow gain would be expected at these frequencies. Figure 
62 shows a peak in the oxidizer side at 2600 Hz which is quite close. On 
the fuel side (Fig. 60) the gain is increasing for increasing frequencies 
below the resonant frequency, which is also expected. The plots of gain 
versus inertance for both the oxidizer and fuel sides (Fig. 61 and 63) are 
also consistent. On the fuel side, increasing inertance lowers the Helmholtz 
frequency and therefore, causes an increase in gain (Fig. 61). On the oxi- 
dizer side, the actual resonant frequency was very near the input frequency 
(2600 Hz). Therefore, either an increase or decrease in inertia would move 
the resonant frequency away from the input frequency causing a reduction in 
gain (Fig. 63). 


To further verify the hypothesis that the Aerojet pie manifolds (especially 
the oxidizer manifold) ■were acting like Helmholtz resonators in resonance 
with the imposed pressure oscillation frequency and, thus, contributing to 
the relatively high gains associated with that injector, a final computer 
run was made in which the volume of the oxidizer pie manifold was decreased 
by a factor of 10. Since this would increase the Helmholtz resonant frequency 
by over a factor of three, a significant reduction in flow gain would be 
expected. The results (see Table X) showed that the gain dropped to 0.12 - 
0.24 of the value with the correct pie manifold volumes, depending on which 


of the four flow summations were used. Therefore, the hypothesis that an 
effective Helmholtz resonance occurs involving the pie manifolds is ana- 

. r . 


lytically predicted. 


The results of this run (see Table X) were also in very good agreement with 
results obtained for the Rocketdyne injector with modified ring grooves and 
orifices. Thus, it was concluded that the differences in gain between the 
Aerojet and Rocketdyne injectors are attributable to (1) the presence of the 
pie manifolds in the Aerojet injector (primary), and (2) differences in the 
ring groove areas and orifice inertances/resistances (secondary). 


Another model result observed in Tables VI through IX is that a higher gain 
results from a spinning IT mode than for a standing IT mode. A spinning IT 
mode chamber pressure input for any ring groove is the same amplitude with 
various phases. The phase of a standing IT mode is either at zero or 180 
degrees, but the amplitude varies with the sine of location of the input. 
Therefore, since the average value of a sine function is 0,64 times the 
maximum value, the average input magnitude for a standing mode is 0.64 
times the average input for a spinning mode (the reference pressure in both 
cases is the pressure at the maximum amplitude) . This difference in flow 
gains for a IT standing mode versus a IT spinning is close to the same ratio. 

As previously mentioned, the objective of the sensitivity analysis was to 
determine generalized injector design criteria to preclude hydraulic coupling. 
It is obvious from the preceding discussion concerning the Rocketdyne and 
Aerojet injectors that the sensitivity analyses were not successful in esta- 
blishing generalized injector design criteria. This is due to the complex 
interactions which may exist between the constituent portions of an injector. 
The analyses did result in a number of observations, however: 

1. Any uncertainty in the linearized orifice resistance due to uncer- 
tainties in the analytical calculation of injector AP should not 
appreciably alter the effect observed by the addition of dams to 
the injector. 
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2. If the ring groove flow calculated by the model is significant com- 
pared to the injector flow, dams in the ring groove could signifi- 
cantly affect the ring groove pressure response, and therefore the 
injector flow oscillations. 

3. For injectors in which complex interactions are absent, the injector 
gain can be simply related to linearized orifice resistance (R) , 
orifice inertance (I), and instability frequency (to) through a 

term called the break frequency (w^) which is simply the linear- 
ized orifice resistance divided by the orifice inertance. These 
simple relationships are: 

gain is proportional to 1/R for ui < to, 

b 

gain is proportional to l/lo for to > to^ 

4. The presence of lai'ge manifold volumes in critical locations may 
introduce complicating mechanisms by acting as Helmholtz resonators 
in resonance with the imposed pressure oscillation. Such action 
can greatly increase the injector gain and nullify completely the 
simple gain relationships presented above. The engine hydraulic 
stability model can be used to show the effect of such large mani- 
fold volumes and suggest design changes to avoid adverse effects. 

Therefore, the engine hydraulic stability model must be utilized in order to 
obtain a quantitative, and in many cases even a qualitative, evaluation of 
geometric injector design options, or operating condition changes, on stability. 

To illustrate the use of the model for such purposes mentioned above, refer- 
ence can be made to the results shown in Tables VI through IX. Figure 59, 
for example, appears to indicate that a reduction in fuel injection tempera- 
ture (for the Aerojet injector) from 230°F to 70°F would decrease the flow 
gain, thus providing less stabilization against the spinning IT mode and 
more stabilization against the resurge mode. Caution must be exercised in 
making such a prognostication, however, since Fig. 59 considers only the 
effect of fuel injection temperature on the injector side of the coupling. 



#' 


4 . ! 


I 

f:: 
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The effect the fuel injection temperature has on the combustion side of the 
coupling is, of course, not predicted by the model, and could be overwhelming. 

Of particular interest to the injector designer are the model predictions 
for the Rocketdyne injector. It is recalled that this technology injector 
had only ^’inular fuel manifold dams. No ring groove dams were ever tested 
on eline_ the fuel or oxidizer sides of the injector. According to Table VIII, 
the fuel side gain* is decreased from 0.00441 (lb/sec)/psi to 0.00367 (lb/sec)/ 
psi when the annular fuel manifold dams arc added. This corresponds to a 
percentage reduction in gain of 16=8 percent. If three symmetrical ring 
groove dams were added to each fuel ring** instead of the three symmetrical 
manifold dams, results from the model show (see Table VIII) that the fuel 
side gain* would be 0.00378 (lb/sec) /psi or a reduction of 14.3 percent 
from the nominal gain with no injector fix. Thus, for the fuel side of the 
Rocketdyne OME technology injector, the addition of manifold dams is predicted 
to be a bit more effective than the addition of ring groove dams in preventing 
coupled instability. However, as shown in Table VIII, simultaneous employment 
of both manifold and ring groove dams would reduce the gain from 0.00441 
(lb/sec)/psi to 0.00266 (lb/sec) /psi . This corresponds to a reduction in gain 
of 39.7 percent. 

The advisability of adding ring groove dams to the oxidizer side of the 
Rocketdyne injector was also explored using the model. These results are 
shown in Table IX. The addition of three symmetrical dams to each ox ring 
(at common theta planes***) is shown in Table IX to reduce the ox side gain 
from 0.00370 (lb/sec>/psi to 0.00346 (Ib/sec) /psi . This 6.5% reduction in 
ox-side gain is less than half the percentage reduction predicted for the 
fuel side. 


* Using the total vector injector flow proportioned by P c amplitude. 

** Dams positioned in each fuel ring at same theta planes as manifold dams 
(see Fig. 35). 

*** See Fig. 38. 
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No effort was made in this program to investigate the effect of ring groove 
dam location on the predicted gain. Such effect may be of considerable im- 
portance*, however, and can be easily studied using the engine hydraulic sta- 
bility computer model. 


* For instance, the location of the ring groove dams in the XRL injector was 
analytically shown to be most critical (Ref. 4). 
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SECTION VI • 


CONCLUSIONS AND RECOMMENDATIONS 

The engine hydraulic stability computer model developed during this program 
is deemed to be of sufficient worth to be included among those models com- 
monly used to investigate combustion instability in rocket engines. The 
range of model applicability is summarized in Table XII. The engine hydrau- 
lic stability model is de.signed to specifically investigate coupling between 
the combustion chamber and the injector hydraulics. It may, however, be input 
in such a fashion to permit the open-loop analysis of feed system hydraulics 
as well. 

The engine hydraulic stability computer model was successful in predicting 
that injector hardware fixes applied to each of three correlation injectors 
would increase combustor stability with respect to the instability mode actu- 
ally observed. 

The model is extremely well suited to: 

1. Quantitatively evaluate the effect of proposed injector fixes on 
stability 

2. Quantitatively evaluate (from a stability standpoint) geometric 
design options or operating condition changes relative to each 
other 

3. Determine potential optimum locations for injector fixes or modi- 
fications such as dams 

4. Determine the depth of penetration of injector face oscillations 
into the feed system or manifold 

The open- loop model is unable to predict stability per se. However, the model 
can be used to look at the injector design detail and indicate what changes 
to make to reduce a high value of injector gaun at expected instability fre- 
quencies . 
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TABLE XII. RANGE OF MODEL APPLICABILITY 


Parameter 

Evaluation 
Range of 
Interest 

Range of 
Correlation 
Test Cases 

Model 

Limitation 

Propellants 

acid/ amine 
LOX /amine 
acid/ 

hydrocarbon 

LOX/ 

hydrocarbon 

IRFNA/UDMH 

N 2 0 4 /MMH 

Liquid propellants 
only 

Thrust 

25 to 50,000 
lbf 

6,000 to 42,000 lbf 

. 

None 

Chamber 

Pressure 

100 to 1000 
psia 

125 to 950 psia 

None 

Mixture 

Ratio 

Maximum per- 
mance ± 20 % 

3.99 (IRFNA/UDMH) 
1.65 (N 2 0 4 /MMH) 

None 

Fuel 

Temperature 

40°F to 0° 
subcooling at 
inj ected 
conditions 

200°F (MMH) 
65 °F (UDMH) 

Liquid propellants 
(no two-phase in- 
jector flow) 

Acid 

Temperature 

40°F to 0° 
subcooling at 
inj ected 
conditions 

65°F (N 2 O 4 ) 
65°F (IRFNA) 

Liquid propellants 
(no two-phase in- 
jector flow) 

LOX 

Temperature 

-298°F to 0° 
subcooling at 
injected 
conditions 


Liquid propellants 
(no two- phase 
injector flow) 

Orifice 

Size 

0.020 to 0.040 
0 .040-in. dia 

0.020- to 0.073-in. dia 

None 

Orifice AP 

25 to 50 psid 
(or 15% P ) 

44 to 184 psid 
(15 to 50% P c ) 

None (no local 
cavitation) 

Frequency 
range of 
coupling 
problem 

100 to 3000 Hz 

1300 to 2600 Hz 

Injector description 
limited to 8 lumps 
per wavelength for 
optimum accuracy 
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It is not recommended that the model be used in an attempt to predict such 
things as a preferred mode or the effect of propellant injection temperature 
since these quantities affect both the combustion side as well as the injector 
hydraulics. The combustion-side response is, of course, input like a boundary 
condition in the model in order to calculate the injector response. The effect 
of variables on the combustion side response is, therefore, not calculated in 
any manner whatsoever and must be assumed. 

The sensitivity analyses were not successful in establishing generalized in- 
jector design criteria to preclude hydraulic coupling. This is because of 
the complex interactions which may exist between the constituent portions of 
the injector. The analyses were most informative, however. For injectors in 
which complex interactions are absent, the injector gain can be simply related 
to orifice resistance, orifice inertance, and instability frequency through 
a term called the break frequency which is simply the orifice resistance 
divided by the orifice inertance. These s^mnle relationships are: 

gain is proportional to. 1/R for w < to, 

b 

gain is proportional to 1/lto for to > to, 

b 

The presence of large manifold volumes in critical locations may introduce 
complicating mechanisms by acting as Helmholtz resonators in resonance with 
the imposed pressure oscillation. Such action can greatly increase the in- 
jector gain and nullify completely the simple gain relationships presented 
above. The engine hydraulic stability model can be used to predict the ef- 
fect of such large manifold volumes and suggest design changes to avoid ad->- 
verse effects. 

The attainment of additional experimental data with which the model can be 
correlated is deemed to be of most urgent concern. In particular, the at- 
tainment of data for which such injector parameters as manifold volume, 
ring groove volume, feed passage inertance, and dam location have been var- 
ied is considered necessary. Such experimental effort would best be conducted 
in conjunction with parametric studies of the experimental variables using 
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the engine hydraulic stability model. Results from the model could then be 
used to help guide the experimental effort. 

Therefore, it is highly recommended that an experimental evaluation of in- 
jector/combustion-coupled instability be conducted* in conjunction with addi- 
tional sensitivity analyses and parameter variation studies. Such effort 
would include the design and fabrication of a versatile injector, testing of 
that injector to determine the most effective means of preventing injector 
hydraulic-coupled instabilities, and attendant analysis and evaluation to 
improve the analytical computer model and develop design criteria and 
recommended procedures. The injector should have sufficient flexibility 
to permit configurational changes (such as the capability for readily in- 
stalling or removing ring groove dams, manifold baffles, orifices in injector 
feed passages, resonators, and the inlet configurations) which are predicted 
to substantially effect the oscillatory flow behavior of the injector but, 
nevertheless, be representative of good injector design practice and fabri- 
cation techniques. 


* Such an effort was originally proposed as Phase II of the Engine Hydraulic 
Stability Contract (NAS9-14801) and is detailed in Ref. 15. 
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APPENDIX A • 


DEVELOPMENT OF GOVERNING EQUATIONS 
FOR LIQUID ROCKET FEED SYSTEMS 


GENERAL FUNDAMENTAL EQUATIONS OF 
CONTINUITY AND MOMENTUM 

The following discussion relates to the development of the differential 
equations governing viscous fluid flow as applied to the analysis of 
liquid rocket feed system stability. Using tensor notation, the contin- 
uity equation describing general fluid motion is (Ref. A-l) : 


+ V ' Cpv) = 0 . 

where 

p is density (mass /volume) 
v is the velocity vector (distance/ time) 
t is time. 


(A-l) 


From isentropic relations, for compression and expansion 



where 

P is pressure (force/area) 

a is the sonic speed of the hydraulic fluid (distance/ time) . 

Thus substituting Eq. (A-2) into Eq. (A-l) will yield a continuity equation 
of 

-^2 + 7 ' ® (A-3) 

a 

Again, using tensor notation, and ignoring body forces, the equation 
governing momentum conservation can be written as (Ref, A-l) 


where 


h .is the momentum head loss vector due to viscous 
surface forces (force/volume) 

Equations (A-3) and (A-4) above are the generalized equations describing 
viscous fluid flow. A few assumptions are introduced to allow the use of 
present mathematical techniques in the solution of the set of equations. 

These assumptions are: 

a. Assume that the system undergoes an initial transient period 
followed by steady state. Thus, only the steady-state solu- 
tion will be sought. 

b. Assume that the steady<-state solution can be characterized 

as being composed of time average values with small perturbed 
values superimposed on top. This is the perturbation assumption. 
Thus let 


p = p + p 

(A-5a) 

X = % + X 

(A- 5b) 

p = p + p 

(A-5c) 

= _h + 

(A-5d) 


where (— ) denotes time averaged values 

( ) denotes small perturbed values which 
approximately equal zero. 

Additionally assume that: 

c. The time averaged velocity vector is also small, 
although 7 is greater than v, and 

d. That the fluid is incompressible, that is the 
time averaged density, p, is constant. 


Assumptions a, b, c, and d above, reduce Equations (A-3) and (A-4) to a 
time averaged steady-state form for continuity 


(A-6) 


V-v. = 0 

and for momentum conservation 

Vp = £ (A— 7 ) 

and, finally, a perturbed steady-state form for continuity 

At 

-J || + p v.£= o (A-8) 

a 

and for momentum conservation 
— 3v *** <*» 

P - = -V P + h (A- 9) 

ot — * 

It should be pointed out that p and a are the only constants in Eqs. (A-6) - 
(A-9) and their values are known for most liquids. 

DEVELOPMENT OF THE VISCOUS 
MOMENTUM HEAD LOSS VECTOR 

Until now no attempt has been made to quantify the viscous momentum head 
loss vector, h. The following discussion will begin this development for 
both general multi-dimensional flow and one-dimensional axial flow in pipes. 


General Multi-dimensional Flow 

It has been shown in Bird (Ref. A-l) that the viscous momentum vector for 
laminar flow can be represented in tensor form as 



(A-10) 


where x is the viscous surface force tensor. 


A-3 


The expressions for t as developed in Schlichting (Ref. A-2) for Newtonian 
fluids are: 


T = 

XX 

-2y 

Sv „ 

ir + t» 

(V-v) 

(A-lla) 

T 

. yy 

-2y 

3y 3 

(V-y) 

(A-llb) 

T = 

ZZ 

-2y 

9v : 

z , 2 

i, 

(V.y) 

(A-llc) 

T = 

xy 

T 

yx 

/ 9v 

" y (ay 

3v \ 

(A-lld) 

T = 

T 

( 3v 

= _y [-Z 

3v \ 
+ . Z 

(A-lle) 

yz 

zy 

\a z 

3y / 

T = 

zx 

T 

xz 

/ 3v 

- - Ur 

3v \ 

(A-llf ) 


where y is the dynamic viscosity (mass/time-distance) 
and x, y, and z denote direction 

Using Eq. (A-10) and assumptions (a), (b) , and (c) , a steady-state time- 
averaged-representation and a steady-state perturbed representation of h 
can be obtained. However, this will not be done here since Eqs. (A-10) and 
(A-ll) would add too great a complexity in the solution of Eqs. (A-3) and 
(A-4) to be treated in this discussion. 

One-dimensional Axial Elow in a Pipe 

For steady-state incompressible isothermal flow in pipes of constant cross- 
sectional area, the Fanning or Darcy equation (Ref. A-3) has shown that the 
pressure changes according to 



D 


AP = 
P 


2 


(A-12) 


where 


| | denotes absolute value 

AP is the pressure change across the length of 
pipe section in the positive axial direction 

f is the fanning friction factor 

D is the hydraulic diameter 

L is the length of pipe section (always positive) 

v is the axial fluid velocity in the positive 
x-direction 

The Fanning equation above has been written in such a fashion as to show 
that in the positive x-direction, the pressure will decrease for positive 
axial fluid velocities, and increase for negative axial fluid velocities. 

Using assumptions (a), (b), and (c) , it can be shown that Eq. (A-4) will 
be reduced to the one-dimensional non-perturbed form of 

3P = . 

3x x (A-13) 

Since these assumptions state that transient terms can be ignored along 
with the non-linear velocity term. Substituting Eq. (A-12) into Eq. (A-13) 
yields an expression for the momentum head loss of 

. ~ 2fp KK 

x D (A-14) 

Now substituting Eqs. (A-5b) , (A-5c) , and (A-5d) from assumption (b) into 

Eq. (A-14) gives for a time averaged head loss term, 

-2fp |v |v 

— _ ' x 1 x 

x D (A-15) 


A-5 



and for a perturbed head loss term, 

r-4fp |v | v 
h = x 1 x 

x D (A-16) 

Finally, using Eq. (A-15), the one-dimensional representation of Eq. (A-7), 
and remembering that AP is opposite in sign from v will show that the 

X 

perturbed head loss term can also be expressed as, 


-2 AP ~ 


v 

T — X 

L v 

(A-17) 


BOUNDARY CONDITIONS FOR 
FUNDAMENTAL EQUATION 

This discussion relates to the boundary equations needed to solve Eqs . (A-6) 
through (A-9) • The spatial flow parameters of pressure and velocity on a 
boundary can be characterized as consisting of the sum of two functions; 
one time dependent, the other time independent. Mathematically, this can 
be represented as 

B Bc = B Bc^ X,y,Z ' ) + B Bc^ X ’ y ’ Z,t ^ (A-18) 

where 

Bg c is the pressure or velocity boundary condition 
at a surface 

B^.fejy.z) is some time Independent function 

DC 

B (x,y,z,t) is some time dependent function 

No boundary conditions are needed for the viscous momentum head loss vector, 
li, since it is a function of velocity (see Eqs. (A-10), (A-ll) , and (A-14) . 

Then from Eq. (A-18) the time averaged boundary conditions used in Eqs. (A-6) 
and (A-7) can be expressed as 

B Bc = B Bc (x,y,z) (A_19) 


A-6 


and the perturbated boundary condition used in Eqs. (A-8) and (A-9) can 
be written as 



= B Bc (x,y,z,t) 


(A-20) 


Finally, since assumption (a) states that the solution is a steady-state 
representation, the initial conditions given for Eqs. (A-8) and (A-9) at 
t = 0 can be taken as equal to zero. Mathematically this is expressed as 


B Ic (t = 0) = 0 


(A-21) 


where 

is the perturbed pressure or velocity initial condition 

The importance of feed-system stability lies in the development of the 
perturbed steady-state equations. Therefore, continued development of the 
general time average steady-state equations, Eq. (A-6) , (A-7) , and (A-19) , 
will be dismissed from further discussion. 


MULTI-DIMENSIONAL INVISCID PERTURBED FLOW 
(THE WAVE EQUATION) 

In the previous sections, a complete multi-dimensional perturbed flow re- 
presentation can be given by Eqs. (A-8) and (A-9), and the perturbated 
description of Eqs. (A-10) and (A-ll) . However, as was stated before, 

Eqs. (A-10) and (A-ll) will add too great a complexity to the solution of 
Eqs. (A-8) and (A-9). Therefore, in the case of general multi-dimensional 
fluid flow, assume that any viscous effects are negligible on the perturbed 
flow, i.e., 

h^fti 0 (A-22) 

Then, using Eq. (A-22), it can be shown that Eqs. (A-8) and (A-9) can be 
combined to yield. 


1_ 

2 

a 


2 

8 P 



2 ^ 

V P 


A-7 


(A-23) 


and that Eq. (A-9) alone reduces to 


p-§f = -VP (A-24) 

Equation (A-23) above is a partial differential equation more commonly known 
as the "wave equation". It can be solved by the method of separation of 
variables, or by the use of Green’s function. Once the solution to Eq.(A-23) 
is known, the time history of the perturbed velocity can be found directly 
from Eq. (A-24). 


ONE-DIMENSIONAL VISCOUS PERTURBED FLOW 

fO 

The perturbed one-dimensional axial flow momentum head loss vector, h , 
given by Eqs. (A-16) or (A-17) is of simple enough form that viscous effects 
can easily be included in one-dimensional flow representations. Expressing 
Eq. (A-8) in one-dimensional form yields for continuity 


«** 3 v 

1 3P . jr X 

2 3t + p 3 x 


= 0 


(A-25) 


By substituting Eq. (A-17) into Eq. (A-9), the one-dimensional momentum 
equation becomes 


Sv x = -3 P 2 [ AP | 

3t 3x ' L j- | V x 
x 


(A-26) 
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APPENDIX B 


CLASSICAL CHAMBER ACOUSTICS 
OF CYLINDRICAL COMBUSTORS 

INTRODUCTION 

When studying liquid rocket feed system dynamics it is necessary to know 
the pressure profile at the injector face. This pressure profile is needed 
for use as a boundary condition on the hydrodynamic feed system differential 
equations. To know how the pressure profile varies across the injector face, 
means that the dynamics of the combustion chamber itself needs to be known. 
This by no means is an easy task since in many cases the dynamics of the 
combustion chamber are coupled with the dynamics of the feed system. But 
assuming that the combustion chamber and feed system are not dynamically 
coupled, a first approximation to the pressure profile on the injector face 
may be found. The discussion which now follows is the simplified or classical 
development of cylindrical combustion chamber acoustics. It will begin by 
first defining a chamber geometry and coordinate system. Once this is done, 
simplifying assumptions for the analytical derivation will be listed with 
the mathematical solution following, 

GEOMETRY 

As stated previously, the geometry of the combustion chamber is cylindrical 
with' the spatial coordinates placed as shown in Fig. B-l, Every location 
within the combustor can be described by its radial direction (r) , its tan- 
gential direction (0), and its axial direction (x) . The injector face is 
located at x = 0, and the radius of the combustor is given by r . 


B-l 



Figure B-l. Combustor Geometry 


ANALYTICAL DEVELOPMENT 

Before the equation governing fluid flow within the combustor is written, 
it is best to first list the assumptions which go into its formulation. 
These assumptions are given as follows: 

a. Assume that the system reaches steady-state, so that only a 
steady-state solution will be sought. 

b. Assume that this steady-state solution is a perturbed solution 
such that it is composed of small oscillatory perturbations of 
pressure, density, velocity, etc., superimposed on top of their 
time-averaged steady-state values. 

c. The time-averaged steady-state velocity is small. 

d. The fluid is considered incompressible. That is the time- 
averaged density is constant. 

e. The fluid flow is inviscid. 

f. The flow is homogeneously non-reacting. That is, the reaction 
kinetics are infinitely fast such that the reaction takes place 
and is completed at one axial, z location. This axial location 
will be considered at the injector face (x = 0). 

g. Gradients in the axial direction are negligible. 


and 


Using assumptions (a) - (g) , the continuity and momentum equations can be 
combined to yield a partial differential equation for the perturbed chamber 
pressure. This equation is developed in Chapter A and given by Eq . (A- 23) 
as 



2 ~ 

8 P 



2 «*» 

a p 


dr 


+ 


1 3P 
r 3t 



where 

P = 


a 

c 


t = 


r , 9 = 


perturbed pressure (force/area) 

sonic speed of the gases in the combustion chamber 
(distance/ time) 

time 

spatial coordinates 


(B->1) 


The solution to Eq. (B-l) is easily obtained when the functional form of the 
dependent variable can be represented as a product solution of the independent 
variables . 


It has been shown experimentally that two types of pressure waves exist in 
the tangential (0) direction of cylindrical combustors. They are standing 
tangential waves and spinning tangential waves. Figure B-2 below graphically 
shows the difference between the two, for a given radial location and one 
period wave (the first tangential) . 

Figure B-2a shows that a standing wave is characterized by spatial oscillatory 
pressure nodes and anti-nodes, much like a vibrating string. However, Fig. B-2b 
clearly shows that spinning waves have no such spatial nodes. All 0 locations 
experience the full amplitude of the tangential wave. It is because of this 
difference that Eq. (B-l) will need to be solved for both these modes. 


THE STANDING TANGENTIAL WAVE 

For the standing tangential wave it will be assumed that the solution to 
Eq. (B-l) can be represented as 

P = R(r) 0(0) T(t) (B-2) 
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Figure B-2. 


The First Tangential Standing and Spinning 
(at a given radial location with amplitude, 


*01 03 



where 


R is some function in r 
0 is some function in 9 
T is some function in t 


Further, Fig. C-2a shows that the functional form of T can be given as 
T = sin cot 

where 

to = angular frequency of oscillations (radians/time) 

Substituting Eqs. (B-2) and (B-3) into Eq . (C-l) will yield 

(B-4a) 
(B-4b) 

where m = separation constant . 


and 


2 d 2 R , dR 

r 77 + r d? 

dr 


+ 


'ui 


2 2 
r - m 


R = 0 


d 2 © u. 2 n n 
r- + m 0 = 0 

a a 4- 


Solving Eqs. (B-4a) and (B-4b) will then yield a solution to Eq. (B-2) of 


P = 


c- J 
1 'm 


ft)* •■’•(?) 


[ c^ cos m0 + sin mej sin mt 


where 

J = Bessel function of the first kind, order m 
m 

Y = Bessel function of the second kind, order m 
m 

c^, Cg. c^ = arbitrary constants 


(B-5) 


The boundary conditions of a standing wave which are needed to solve for the con- 
stants in Eq. (B-5) are 
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P(r = 0) = finite 


sp 

3r 


= 0 


(B-6) 

(B-7) 


r = r 


w 


and 


P ( 0 =0) = P (6 = 2 it) 


ap 

30 


0 = o 


3P 

30 


■ 0 


(B-8) 


(B-9) 


0 = 2tt 


Applying Eqs. (B-6) and (B-9) to Eq. (B-5) shows that the constants c 2 and 
equal zero. For Eq. (B-8) to hold, the constant "m" must be an integer, so 
that 


m = 0, 1, 2, 3, 


(B-10) 


Finally, applying Eq. (B-7) to Eq. (B-5) will show that the system must oscil- 
late at discrete frequencies such that 


d J 

m 




dr 


= 0 


(B-lla) 


r = r 


w 


where 


n = 1, 2, 3, . . . . 


(B-llb) 


Table B-l below shows some of the solutions to Eq. (B-ll) , giving the eigen- 
values of a) for the various transverse modes. Fig. B-3 then graphs out the 
m,n 

Bessel function, J m , versus the non-dimensional radius, (r/r w ), for the first 
and second tangentials, and the first radial acoustic modes. 


I 
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TABLE B-l. TRANSVERSE ACOUSTIC MODES 


)( 

|jj 

I 


m 

n 

/ w r \ 

1 m,n w \ 



Transverse Mode 

V a c ) 

1 

i 

1.8413 

First tangential 

2 

i 

3.0543 

Second tangential 

0 

2 

3.8317 

First radial 

3 

1 

4.2012 

Third tangential 

0 

3 

7.0156 

Second radial 

1 

2 

5.3313 

Combined first tangential and 
first radial 

1 

3 

8.5263 

Combined first tangential and 
second radial 

2 

2 

6.7060 

Combined second tangential and 
first radial 


Thus, in light of the previous discussion, Eq . (B-5) finally reduces to 



(B-12) 


Where P ^ = c i c 3 (the amplitude of the first radial mode at r = 0) . 


Equation (B-12) shows the complete functional form of all the oscillatory modes 
associated with standing tangential waves. The actual oscillatory pressure for 
a given "r", "0", and "t" is made up of the sum of the contributing pressures 
of each mode so that 


P 



oo 



n-1 


(B-13) 
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Figure B-3. Bessel Functions for the First Three Acoustic 
Modes Plotted Versus Non-Dimensional Radius 
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THE SPINNING TANGENTIAL WAVE 

For the spinning tangential wave it will be assumed that the solution to 
Eq. (B-l) can be represented as 


P = R(r) F(e,t) 

where F is some function in 0 and t. 


Further, Fig. B-2b shows that the functional form of F is given by 
F = sin (u:t + m0) 

Substituting Eqs. (B-14) and (B-15) into Eq; (B-l) will yield 

2 d 2 R , dr , | u 2 2 2 „ 

r — r— + r - 7 - + { — - r - m R = 0. 

,2 dr 2 

dr a 

c 


(B-15) 


(B-4a) 


This is the same result for the radially dependent function that was obtained 
for the standing wave. Thus the solution to Eq. (B-14) is simply 

/V I MY“ MY- I 

P = c, J (- c 0 Y — I sin (u)t + m0) (B-16) 

1 m a 2 . m a I 

The boundary conditions of a spinning wave are the same as those for standing 
waves (Eqs. B-6) - (B-8)), with the exception being that Eq. (B~9) is given 
only as 


i. B-6) 

- (B-8)), with the exception being that Eq. 

(B-9) is given 





r.i 

3P 

SP 



30 e - 

0 = 30 0 = 2TT 

(B-17) 

( 


Therefore once again c ^ equals zero, "m’’ is given by Eq. (B-10) , and the 
system will oscillate at the same frequencies given by Eq. (B-ll) . 


Thus for a cylindrical combustor with spinning tangential waves the 
oscillatory pressure is given by 
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p = p J 
m,n rad m 


( u r \ 

-sH * 1 


in (o) t + m0) 
m,n 


(B-18) 


where here P ^ = c^ (again the amplitude of the first 
radial mode at r = 0) . 


As is easily seen, the only difference between Eq, (B-12) for standing 
waves and Eq. (B-18) for spinning waves is in replacing l ‘cos m0 sin cot 1 ’ 
by "sin (cot + m0)". Therefore Table B~1 and Fig. B-3, which were generated 
from the radial function only, will still apply to both standing as well 
as spinning waves. 
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APPENDIX C 

COMPUTER MODEL 
DOCUMENTATION OF AEROJET 
OME TECHNOLOGY INJECTOR 
FUEL SYSTEM FREQUENCY RESPONSE 
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Figure C-l. Data Deck for Aerojet Injector Fuel Side 
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Angle and Radius for Each Chamber Pressure Input 


Figure C-l, (Continued) 













270 . 1.212 0 . 1.212 120 . 1.212 

240 . .688 0 . .688 120 . .688 

240 . .162 0 . .162 120 . .162 

0 . 0 "' 

0.0 

C 01 R= 16*0. , 15*23. 3, 1 5*0 . , 6* 2 . 02 , 2 1 *0. ,6*. ?? 3, 6*. 287 i6* . 3 54 » 6*. 5 , 

6 *. 761 , 6 * 1 . 33 , 6 * 7 . ,6 *4 1.7 , 0 . 0 , 3 * 0 . 0 , 0 . 0 , 3 * 0 . 0 , 0 . 0 , 3 * 0 . 0 , 

0 . 0 , ?* 0 ., 0 . 0 , 2 * 0 ., 0 . 0 , 2 * 0 . , 0 . 0 , 2 * 0 . , 0 . 0 , 2 * 0 . , 

0 . 0 , 2 * 0 . , 0 . 0 , 0 . , 0 . 0 , 0 . , n . 0 , 0 ., 0 . 0 , 0 . , 0 . 0 , 0 ., • 

' 0 . O', 0 . 0 , 3 * 0 . 0 , 6 * 0 . 0, 12 + 1294 . , 9 * 859 . , 9 * 862 . ,6 *534 . , 6 * 643 . , 

3 * 474 . , 3 * 883 . , 3*24 58 ., 

Z = 16 +. C 046 04 , 15 *. 044 07 , 1 5 * . 03732 , 

6*.0 06 545, .008 75 7, .005 913 , .002 84 1 , .005 69 3, .0 0923 1 , .0 05693 , .009231 , 
. 00075 7 ,. 00591 3 , . 002841 , . 0 0569 3 , . 00923 1 , . 005 6 ^ 3 , . 009231 , . 003757 , 
. 005913, .002 84 1 ,. 00569 3 ,. 009 231 , . 00569 3 C 09 23 1 , 6* .006959 , 6 * . 005274 , 

6 *. 0067 8 8 , 6 *. 0 09 161 , 6 *. 013 42 , 6 *. 02 23 3 , 6 * . 064 8 1 , 6 * . 2637 , 1 2 *. 1 742 , 

~ 9 *. 125 5 , 9 *. 120 2 , 6 *. 173 5 , 6 *. 173 2 , 3 *. 344 7 , 3 *. 2 1 , 3 *. 075 3 1 , 12 *. 03307 , 

9 *. 01 4 29 , 9 *. 01 68 4 , 6 * . 01 36 7 , 6 *. 0165 5 , 3 * . 012 09 , 3*. 02 245 , 3 *. 062 8 6 , 

V = 15 * 2 . 756 ,. IS 13 , 2 *. 23 33 , 3 *. 19 13,?*. 23 S 3 , 3 *. 13 13 , 2 *. 2383 , 2 *. 1613 , 1 . 195 , 
. 5272 , . 1576 , 1 . 195 , . 5272 , . 1576 , 1 . 195 , . 5272 , . 1 576 , 1 . 19 5 , . 5272 , . 1576 , 
1 . 195 , . 5272 , . 1576 , 1 . 195 ,. 5 272 , . 1576 , . 06416 ,. 0642 ,. 0642 , . 06416 , 

.06416 , 2 *. 0642 , 2 *. 0641 6 , 2 *. 064 ?, . 06416 , 2 *. 1192 , . 1192 , 2 *. 1 192 ,. 1192 , 
" 2 *. 11 92 , . 1 19 2 , 2 *. 0394 , . 0994 , 2 * . 0894 , . 0 89 4 , 2 * . 0894 , . 08 94 , . 09369 , 

. 0937, .093 69 ,. 0937 ,. 09 369 , . 0937 ,. 05724 , .0572 , . 05724, .0572 ,. 05724 , 

. 057 ?, 3 *. 1172 , 3 *. 064 14 , 3 *. 0231 8 , 


C=99*47880 , SEND 



Figure C-l. (Continued) 


C-4 




REAL INPUT MATRIX AMPLITUDES - PSI 

0.0 7.071 IE- 0 1 9.6593F-01 9.6593E-01 7.0711E-01 

— 2. 58 8?E-0L -7.071 IE-01 -9.6593E-01 -9.6593E-01 -7.071 IE-01 


2. 5882 F- 01 
-9.2173E-01 
6.9644E-01 
1. 5787E-01 
0.0 
0.0 

-1.3434F- 01 


7.51'iOF-Ol 9. 0 089E- 01 7.5140F-01 1.7033F-01 

-9.2 173E-0 1 >A.o0AAE-0r'"l'.7033E-01*‘ 6. 9644E-01 
1.5787E-01 - 4. 54 57 E— 01 -8.5432E-01 -8.5432E-01 
6.8675F-01 6. 8675E- 01 0.0 -6.8675E-01 

5 . 5602.E- 0 1 5. 56C2F-01 0. 0 -5.5602E-01 

.4. 6390E- 0 1 -2.3195E-01 -2.3195E-01 2.6867E-0t 

6. 3622E- 02 -3.1811E-02 -3.1811E-02 


2.508 2E— 01 
-2.5882E-01 
-4 . 9044b— 01 
9. 091 4 E- 01 
-4.5457E-01 
—6 .8675 E— 0 1 
-5. 5602E-01 
— 1.3434E-01 


I MAGINAP.V 

0.0 

9.6593F 
-9.6593E 
3.3548E 
5.8439E 
— 8.95 3 3E 
-7.9299F 
-6.4204E 
-2.3268E 


INPUT MATRIX 
-7.0711E- 
-01 7.0711E- 

-01 -6. 3050E- 
-01 -3.3548E 
-01 8.95336 

-01 -3. 9650E- 
-01 -3.2102E- 
01 0.0 
-01 0.0 


AMPLITUDES - 
01 -2.5882E- 
2.5882E- 
0.0 

8. 49 47 E- 
7. 8734E- 
3.96 5 0F- 
3.2102E- 
4.0175E- 
5. 5098E- 


01 

01 

■01 

■01 

01 

01 


PSI 

01 2. 5882F-01 

01 -2.5882E-01 
6. 3050E-01 
•01 -9.6599E-01 
01 3.1 095E-01 

01 ’7. 9299E-01 

01 6.42 04 E-0 1 
■01 -4. 0175E— 01 

02 — 5.5098E— 02 


7.07 11 E-0 1 
-7. 07 1 1 E— 0 1 
9. 6599 F-0 1 
-5.8439E-01 
—3 . 109 5 E-0 1 
3. 9650E-01 
3.2102E-01 
0.0 


9. 6593 F-0 l 
—9 .6593 E— 0 1 
8.4947E-01 
0.0 

-7.8 73 4F— 0 1 
— 3. 9650 E— 01 
—3.2102 E-01 
2. 3268E-01 


Figure C-2. 


Real and Imaginary Input Matrix Amplitudes 
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Figure C-3. Input and Output Data for Flowrates 
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Figure C-4. Input and Output Data for Pressure Nodes 
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Figure C-4. (Continued) 
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Figure D-l, Data Deck for Aerojet 
Injector Oxidizer Side 
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Figure D-l. (Continued) 
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Figure D-2, Real and Imaginary Input Matrix Amplitudes 
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Figure D-4. Input and Output Data for Pressure Nodes 


D-9 


42 

1. 170E-01 

3.96 IE+04 

7. 578E-01 

204.9 

42 

43 

94 

95 

142 

43 

1. 170E-01 

3.961E+04 

8. 28CE-01 

347.9 

95 



87 

143 

44 

8.320E-0 2 

3.961E+04 

7.282E-01 

10. 0 

96 



97 

144 

45 

8. 320E-02 

3.96 IE+04 

7. C05E-01 

32 3. 8 

44 

45 

97 

98 

145 

46 

8. 320E-02 

3 . 96 1 F +0 4 

6. 89 5E-0 l 

112.5 

98 



99 

146 

47 

8.320E-02 

3. 96 IE 4-0 4 

7. 334E-01 

129.9 

99 



100 

147 

48 

8. 320 E- 02 

3.961E+04 

7. 155E-01 

85. 7 

47 

100 

46 

10 1 

148 

49 

8. 320E-U2 

3 . 96 IE 404 

fc. 856E-C1 

232.4 

101 



102 

149 

5Q 

8, 320E-02 

3.96 IE 40 4 

7.282E-01 

251. 0 

102 



103 

150 

51 

8. 32 0 E- 0 2 

3.961E4-04 

6. 584F.-01 

204. 6 

48 

49 

103 

10 4 

151 

52 

8. 320 E- 02 

3.96 IE 4-04 

6.985E-01 

352.3 

104 



96 

152 

53 

1.430E-01 

3. °6 IE+04 

5. 895E-0 1 

327. 2 

50 

105 


106 

153 

54 

L.480E-01 

3.96 IE *04 

4. 180E-01 

303.4 

51 

l C6 


10 7 

154 

55 

1.4 8 OF— 0 1 

3 .96 IE +04 

6. 054E--01 

88. 7 

52 

107 


108 

155 

56 

1.480E-01 

3 .96 IE +0 4 

4. 313F-01 

7 1.0 

53 

ICP 


109 

1 56 

57 

1.48 0E-0 1 

3.961E+04 

5. 529E-01 

20 8. 7 

54 

109 


no 

157 

58 

1.480E-01 

3 .96 IE +04 

3. 815E-01 

189.1 

55 

no 


105 

158 

59 

1.01 0E-01 

3.96 IE +0 4 

3.16CE-01 

33 5. 4 

56 

1 11 


112 

159 

60 

1.010E-01 

3 .96 IE +0 4 

2. 027E-01 

310.4 

57 

112 


113 

160 

61 

l.OlOE-O'l 

3.961E+04 

3.43 ?E-0 1 

97. 6 

58 

113 


114 

161 

62 

1.010E-01 

3. 96 IE + 04 

l. 938E-01 

74. 5 

59 

114 


115 

162 

63 

i.OlOE-Ol 

3.96 IE +0 4 

3. 479E-01 

208. 8 

60 

115 


116 

163 

64 

1.010F-01 

3.96 IE+04 

6. 123E-01 

176. 9 

61 

116 


111 

164 

65 

1.230E-01 

3.96 IE +0 4 

4. 44 8E-01 

326. 7 

62 

63 

L 17 

118 

165 

66 

1.230E-01 

3.961E+04 

2. 60CE-01 

82.4 

64 

65 

118 

119 

166 

67 

1.230E-01 

3.961E+04 

2.5 5 IE- 01 

198. 7 

66 

67 

119 

117 

167 

68 

8.020E-02 

3. 96 IE +04 

2. 379E-01 

192. 6 

68 

69 

120 

121 

i6e 

69 

8.020E-02 

3.96IE+04 

2. 40 IF— 01 

7 8.3 

70 

71 

121 

122 

169 

70 

8. 020 E- 02 

3 .96 IE +0 4 

2.432E-01 

316.5 

72 

73 

122 

120 

170 

7 1 

1. 520E-02 

3.96 IE+04 

5. 520E-03 

292.9 

74 



17 1 



AEROJET CME INJECTOR 
OXIDIZER SIDE, NO RING DAMS 
SPINNING FIRST TANGENTIAL MODE 


Figure D-4. (Continued) 


APPENDIX E 

COMPUTER MODEL DOCUMENTATION OF ROCKETDYNE 
OME TECHNOLOGY INJECTOR 

FUEL SYSTEM FREQUENCY RESPONSE 


E-l 


ROCK ET DYNE OME INJECTOR 
FUEL SIDE, NO RING DAMS 
STANDING FIRST TANGENTIAL 


78 

39 

12 


1 

130 

-10 

-28 

2 

131 

- 1 1 

-29 

3 

132 

-1 2 

-30 

4 

133 

-13 

-31 

5 

134 

-14 

-32 

6 

135 

-15 

-33 

7 

136 

-16 

-34 

a 

137 

-17 

-35 

9 

138 

-18 

-36 

10 

-19 



u 

-20 



12 

-21 



13 

-22 



14 

-23 



15 

-24 



16 

-25 



17 

-26 



18 

-27 



19 

100 

-101 

-139 

20 

101 

-102 

-140 

21 

102 

-103 

-141 

22 

103 

-104 

-142 

23 

104 

-105 

-143 

24 

105 

-106 

-144 

25 

106 

-107 

-145 

26 

107 

-108 

-146 

27 

108 

-100 

-147 

28 

-37 



29 

-38 



30 

-39 



31 

-40 



32 

-41 



33 

-42 



34 

-43 



35 

-44 



36 

-45 



37 

109 

-110 

-148 

38 

110 

-11 l 

-149 

39 

111 

-112 

-150 

40 

112 

-113 

-151 

41 

113 

-114 

-152 

42 

114 

-115 

-153 

43 

115 

-116 

- 154 

44 

116 

-117 

-155 


MODE 


-46 

- 64 - 

-73 

-131 


-4 7 

-48 

-65 

-74 

-132 

-49 

-66 

-75 

-133 


-50 

-67 

-76 

-134 


-51 

-52 

-68 

-77 

-135 

-53 

-69 

-78 

-136 


-54 

-7 0 

-79 

-137 


-55 

-56 

-71 

-80 

-138 

-57 

-72 

-81 

-130 



Figure E-l. Data Deck for Rocketdyne OME 
Technology Injector Fuel Side 
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Figure E-l. (Continued) 
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Figure E~l. (Continued) 
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Figure E-2 , Real and Imaginary Input Matrix Amplitudes 
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Figure E-3, Input and Output Data for Flowrates 
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P ( 47 ) 

1 .06CE *01 

9. 060E-02 

3 .7761-05 

94.7 

49 

P( 

3) 

P ( 47 ) 

5. 2S0E+00 

4.530F-02 

. 1.932F— 04 

98.4 

50 

P( 

4) 

P ( 48 ) 

5.250E+00 

4. 530F-02 

1.287E-04 

278.7 

51 

P( 

5) 

P ( 43 ) 

1 . 06CE *01 

9.060E-02 

1.367E-05 

111. 1 

52 

P 1 

5) 

P ( 49 1 

1.060E+01 

9.060E-02 

1.423E-04 

278.6 

53 

P( 

6) 

P( 49) 

5.2 90E *00 

4. 530E-02 

2.306E-04 

278.0 

54 

P( 

7) 

P ( 50 ) 

5.2S0E+00 

4. 530E-02 

1.188E-04 

277. 5 

55 

P( 

8) 

P( 50) 

1. 0 6CE *01, 

9. C60E-02 

l . 3 38F-04 

279. 5 

56 

P( 

8) 

P ( 5 l ) 

1 .060E *01 

9.060E-02 

1.046F-04 

100. 1 

57 

P ( 

9) 

P ( 5 1 ) 

5.25CE+C0 

4.530E-02 

3. 741E-05 

95.0 

58 

P ( 46 ) 

P ( 52 ) 

3.52CF *00 

3. 030E-02 

2.772F-04 

97.8 

59 

P (47 ) 

P( 53) 

3 • 520F + 00 

3. 030E-02 

l. 68UF-04 

97.4 

60 

P ( 48 ) 

. P ( 54 ) 

3.52CE+00 

3. 030E-02 

7.992F-05 
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61 
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P( 55) 

3.520E*00 

3. 030E-02 
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3.52CE+00 

3.030E-02 
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278. 3 
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3.030F-02 
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1.231E— 05 
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PRESSURE 

: VOLUME 

AC. VEL. 

MAGNITUDE 

NODE 

CU IN 

IN/SEC 

PSI/PSI 

I 

1.34CE+00 

5 . 040E *04 

2. 057E-01 

2 

1.34OE+0C 

5.040E+04 

2. 106E-01 

3 

1 . 34 CE + 00 

5.04CE+04 

1.232E-01 

4 

1.340E+0C 

5 • 040F + 0 4 

5.65 BE- 02 

5 

1.340E+00 

5.040E+04 

1. 716E-0L 

6 

1.340E+00 

5.040E+04 

2.1 16E-0 l 

T 

1. 340E+0C 

5.040E+04 

l. 492E-01 

8 

1.340E+00 

5.040E+04 

2.891E-02 

9 

1.340E+00 

5.040E+04 

8. 836E-02 

10 

4.430E-02 

5. C40F *-04 

3. 81CF-01 

11 

4.430E— 02 

5.040E+04 

3. 854E-01 

12 

4.43CE-02 

5.040F+04 

2. 137E-CI 

13 

4.430E-02 

5.040F*04 

8. 527E-02 

14 

4.430E-02 

5.04CF+04 

3. 142E-01 

15 

4 . 43 0E- 02 

5. 040E+04 

3.995E-01 

L6 

4. 430E-02 

5.040F+04 

2.957E-CI 

17 

4.43 CE-02 

5.040E+04 

7. 122E-02 

18 

4.43CE-02 

5. 04GE+04 

1.85 EE -01 

19 

4. 340E-02 

5. 040E + 04 

5. 355E-CI 

20 

4.340E-02 

5.040F+04 

5.393E-01 

2 1 

4.340E-02 

5.040E+04 

2.926F-01 

22 

4.343E-02 

5.040E+04 

1. 093E-01 

23 

4. 340E-02 

5.040E+04 

4 . 397E-0 1 

24 

4.340E-02 

5.040E+04 

5.657E-C1 

25 

4.340E-02 

5.040EF04 

4.262E-01 

26 

4. 340 E- 02 

5.040E«-04 

9. 966E-02 

27 

4 . 340E-0 2 

5.040E+04 

2.732E-01 
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S.410E-02 

5.G40E+04 

3. 67 CE-01 

29 
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5 . 040E+04 

3.716E-01 
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5.040E+04 
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Figure E-4. Input and Output Data for Pressure Nodes 
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Figure G-3, Input and Output Data for Flowrates 


G-5 


48 

P( 121 

PIN! 4) 

2. 030E+C2 

2.420E-02 


1.459E-03 

279.9 

49 

P( 131 

PIN! 5) 

2.0^0E+02 

2.420E-02 


2.251E-03 

300.3 

50 

P( 14 1 

PIN! 6) 

2. 0306 +02 

2. 420E-02 


2. 5056-03 

i24. 3 

51 

PI 15) 

PIN! 7) 

2.030F+02 

2. 42GE-02 


1.482E-03 

IC3.0 

52 

P ( 16) 

PIN! 8) 

2.0 30E + 02 

2.420E-02 


1.308E-03 

100. 0 

53 

P< 17) 

PIN! 9) 

2.030E+C2 

2. 42 OE -02 


2.1 50E-03 

120. 3 

54 

PI 18) 

P INI 10) 

2. 030E+02 

2 .420E-02 


2.569E-03 

20 3. 6 

55 

PI 19) 

PIN! 11) 

2. 030E+02 

2. 420F-02 


1.544E-03 

282. 4 

56 

PI 20) 

PIN! 12) 

2 . 030E *-02 

2.420F-02 


1.459F-03 

279.9 

57 

P 1 2 1 ) 

P INI 13) 

2. Q30E+02 

2. 420E-C2 


2.251 E-03 

300. 3 

58 

P I 22 ) 

PIN! 14) 

2 • C3CE+02 

2.420E-02 


2. 505 E- 03 

124.3 

59 

P (23) 

PIN! 15) 

2.02GE+G2 

2.420E-02 


1.482E-03 

103.0 

60 

PI 24) 

P INI 16) 

2. 02CE+02 

2. 420E-02 


1 .386E-03 

1C0.0 

61 

P I 25 ) 

P INI 17) 

2.2 10E + 02 

2. 600E-02 


2.217E-03 

123. 6 

62 

PI 26) 

PIN! 18) 

2.2 10E + 02 

2. 6 COE -02 


2.216E-03 

300.2 

63 

PI 27) 

PIN! 19) 

2.2 10E+02 

2 • 600E-02 


1.664E-03 

287.5 

64 

PI 28) 

P IN 1 20 ) 

2.2106+02 

2. 6 OOF- 02 


2.27 IE-03 

203.4 

65 

PI 29) 

P INI 21) 

2.2 10E + 02 

2. 6CCE-02 


2. 182E-03 

120. 7 

66 

PI 30) 

P IN I 22) 

2.2 10E+02 

2. 600E-02 


1.627 E-03 

10 7.8 

67 

P I 3 1 ) 

PIN! 23) 

2.2 10E + 02 

2. 6CCE-02 


2. 21 7 E-03 

123.6 

68 

PI 32 1 

P IN I 24 ) 

2.2 10E+02 

2. 600E-02 


2.2 16E-0 3 

300. 2 

69 

P (33) 

P IN (25) 

2.2 10E + C2 

2. 600E-02 


1.664E-03 

287. 5 

70 

P< 34) 

P IN I 26) 

2.2 10F+02 

2.600E-02 


2. 271 F- 03 

303.4 

71 

P (35 I 

PIN! 27) 

2.2 10E+02 

2. 600E-02 


2.182E-03 

120.7 

72 

P I 36 ) 

PIN! 28) 

2.2 10E+02 

2 • 6C0E-02 


1.627E-03 

107. 0 

73 

PI 5) 

PI 37) 

9.D00E-02 

6.5006-03 


l . 7 53E-03 

106.3 

74 

PI 6) 

P I 38 ) 

9. 0CGF-02 

6. 500E-03 


1. 496E-03 

289. 6 

75 

PI 7) 

PI 39) 

9. OC3F-02 

6.500E-03 


1.753E-03 

106.3 

76 

PI 8) 

P 1 40 > 

9.000E-02 

6. 5C0E-03 


1.496E-03 

289. 6 

77 

P I 38 ) 

P I 4 1 ) 

0.0 

1.756F-01 


3.446E-04 

114. 1 

78 

P I 4 1 ) 

P I 42 ) 

0.0 

i. 756E-01 


5.1 87E-05 

294. 1 

79 

P 1 42 ) 

P (43) 

0.0 

1.756F-01 


7.808E-06 

114.1 

80 

P (43) 

PI 44) 

0.0 

1.75 6E-01 


1. 175E-06 

294. 1 

81 

P I 44 ) 

P (45 ) 

0.0 

1. 756E-01 


1.768E-07 

114. 1 

82 

P I 45 ) 

P I 46 ) 

0.0 

1.756E-01 


2.bC2E-08 

294.1 

83 

P 140 ) 

P I 47 ) 

0.0 

1 • 756E-01 


3.446E-04 

114. 1 

84 

P I 47 ) 

PI48) 

0.0 

1. 756E-01 


5. 187E-05 

294 . 1 

85 

P I 48 ) 

P I 49 ) 

0.0 

1. 756E-01 


7.808E— 06 

114. 1 

86 

P 1 49 ) 

P I 50 ) 

0.0 

1. 756E-01 


1. 175E-06 

294. 1 

87 

PI 50) 

P I 5 l ) 

0.0 

1. 756E-01 


1.76BE-07 

114.1 

88 

P I 5 1 ) 

P I 52 ) 

0.0 

1. 756E-01 


2.602E-08 

294. 1 

XRL 

INJFCTOR 







FUGL 

SIDE, NO RING 

DAMS 






ANOMOLOUS SECOND TANGENTIAL MOOE 










LB/SEC/PSI 

% FLOW/?; PC 

TOTAL 1 

INJECTOR FLOW* 




5 

.5052E- 02 

1 . 34 1 OE +00 

TOTAL VECTOR INJECTOR 

FLOW= 



5 

.4374C-02 

1. 32456+00 

TOTAL 

INJFCTOR FLOW PROPORTIONED 

BY PC AMPL I TUDES= 

5 

.5052F-0 2 

1 .34106 + 00 

TOTAL VECTOR INJECTOR 

FLOW PROPORTIONED BY PC 

AMPL I TUDES= 

5 

.43 74 E- 02 

1.3245E+00 


Figure G-3 


(Continued) 


PRESSURE 

: VOLUME 

AC. VEL . 

MAGN I TUDE 

PHASE 

FLOSS IN 

FLOWS OUT 

NODE 

CU IN 

IN/SFC 

PSI /PSI 

DEGRFES 





l 

3. 200E+00 

4.770F+04 

l. 500E-01 

23.7 

l 


5 


2 

3. 200E+0C 

4.770F+04 

1. 484E-C1 

204.2 

2 


6 


3 

3.200E+0C 

4. 770E + 04 

1. 500E-01 

23.7 

3 


7 


4 

3.200E+00 

4. 770E t-04 

1.484E— 01 

204.2 

4 


0 


5 

3.025E+00 

4.77GE+04 

3.640E-01 

23.3 

5 


13 

73 

6 

3. 025E+OC 

4.770F+04 

3. 505E-01 

205.2 

6 


14 

74 

7 

3.025E+0C 

4. 770F+04 

3.640E-01 

23.3 

7 


15 

75 

8 

3.025E+0C 

4 .7 73E+0 4 

3. 505E-01 

205.2 

8 


16 

76 

9 

l . 0 6 OE +0 0 

4 . 770E+04 

4.436E-01 

21.5 

9 

17 

18 

45 

10 

7.540E-01 

4. 770E+04 

3.270E-01 

208. 1 

1 8 


19 

46 

11 

7. 54 OE— 0 1 

4.770F+04 

6.78 7E-01 

200.8 

19 


20 

47 

12 

7.54CE-01 

4 .770E+04 

7.078E-01 

200.0 

20 


21 

4 8 

13 

1.060E+00 

4.770E+04 

4. 215E-01 

203.9 

21 

10 

22 

49 

14 

7.540E-01 

4. 770E+04 

3. 357E-01 

25. 0 

22 


23 

50 

15 

7.540E-01 

4.770E+04 

6.C56E-01 

19.4 

23 


24 

51 

16 

7.540E-01 

4.770E+04 

7. 187E-0L 

18.6 

24 


25 

52 

17 

1.060E+00 

4.770E+04 

4.43CE-01 

21.5 

25 

11 

26 

53 

18 

7. 540 F- 01 

4 . 770E+04 

3.27CE-01 

208. 1 

26 


27 

54 

19 

7.54CE-01 

4 .770F. *-04 

6.787E-01 

2C0.8 

27 


2 8 

55 

20 

7.540E-01 

4.770E+04 

7.07CE— 01 

200. 0 

28 


29 

56 

21 

1.060E+00 

4. 770E + 04 

4.215E-01 

203.9 

29 

12 

30 

57 

22 

7. 540E-01 

4.770E+04 

3.357E-C1 

25. 0 

30 


31 

58 

23 

7.540E— 01 

4.77CE+04 

6.856E-01 

19.4 

31 


32 

59 

24 

7.540E-01 

4.770E+04 

7. 187E-01 

18.6 

32 



17 

60 

25 

7. 620 E- 01 

4 . 770E + 04 

3.67 EC- 0 1 

23. 7 

33 

13 

34 

61 

26 

5.060E-01 

4.770E+04 

3.936E-01 

208.2 

34 


35 

62 

27 

5. 060 E- 01 

4 . 7 70E+04 

6.C41E-01 

203.9 

35 


36 

63 

28 

7.620E-01 

4.770E+04 

3.560E-01 

205. 6 

36 

14 

37 

64 

29 

5.060F-0 1 

4.770E+04 

3.976E-01 

26. 7 

37 


38 

65 

30 

5.060E-01 

4. 770E+04 

6. C8SE-01 

22.9 

38 


39 

66 

31 

7.620E-01 

4 • 770F + 04 

3.67 8E-01 

23. 7 

39 

15 

40 

67 

32 

5.060E-01 

4 .770E *0 4 

3. 93EE-01 

208.2 

40 


41 

68 

33 

5.060E-01 

4. 770E + 04 

6.04 IE- 01 

203.9 

41 


42 

69 

34 

7.620E-0 1 

4. 770E+04 

3. 560E-01 

205.6 

42 

16 

43 

70 

35 

5. 060 E- 01 

4.770E»-04 

3.976F-01 

26.7 

43 


44 

71 

36 

5.060E-01 

4.770E+04 

6.08 9E-01 

22. 9 

44 


33 

72 

37 

3.740E*0C 

4.770E«-04 

4.565E-01 

21.8 

73 


9 


38 

3.740E *00 

4.770E+04 

4. 296E-01 

204. 1 

74 


10 

77 

39 

3.740E«-00 

4 . 7 70F *04 

4.56 5E-01 

21. 8 

75 


11 


40 

3. 740E+0C 

4.770E+04 

4.296E-01 

204. 1 

76 


12 

83 

41 

4. 420E+00 

4. 770E + 04 

6. 467E-02 

24. 1 

77 


78 



Figure G-4. Input and Output Data for Pressure Nodes 
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Figure H-l, Data Deck for XRL Oxidizer Side 
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Figure H-l, (Continued) 
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l.OOOOE+OO 

-l.OOOOE+OO 

-l.OOOOE+OO 

IMAGINARY INPUT MATRIX AMPL 
0.0 0.0 

ITUnES - PSI 
0.0 

0. 0 

0.0 

0.0 

0.0 

0. 0 

c.o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0. 0 

0.0 

0.0 

0.0 

0.0 

0 .0 

0.0 

0.0 

0. 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0. 0 

0.0 


INPUT FREQUENCY= 1300.0 


Figure H-2. Real and Imaginary Input Matrix Amplitudes 
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FLOW 

UPSTREAM 

COWNSTREAM 

RESISTANCE 

I NERTANCF 

AMPLITUDE 

PHASE 


PRESSURE 

PRESSURE 

SEC /IN SQ 

SEC SO/IN SQ 

L8/SEC/PSI 

CFGRFE! 

1 

P ( 6 1 ) 

P( 1) 

1 .880F-02 

4.400E— 03 

7.939F-03 

296.3 

2 

P(6l » 

P( 2) 

1.8ECE-02 

4. 400C-03 

7.9 39E-0 3 

116.3 

3 

P( l ) 

P< 3) 

1 .80OE- 02 

4.4C0F-03 

6. 649 F— 03 

296.3 

4 

P( 2) 

P( 4) 

1.88 or- 0? 

4.4 COE — 03 

6.649F-03 

116.3 

5 

P( 3) 

P( 5) 

i.aeoE-02 

4. 4 COE — 03 

4.278E-03 

296.3 

6 

P< 41 

P( 6) 

1 .880F-02 

4. 400E-03 

4.278F-03 

116. 3 

7 

P( 5) 

P( 7) 

1.8 60E-02 

4.400E-03 

1.2 12F.-03 

296.4 

8 

P( 6) 

P< 8) 

1.880E-02 

4.400E-03 

1.212E-03 

116.4 

9 

P( 7) 

P( 9) 

1.880F-C2 

4.400E-03 

2.05 IE-03 

116.2 

10 

P( 8) 

P ( 10) 

1 .880F-02 

4.4C0F-03 

2.05) F-03 

296.2 

11 

P< 9) 

P( 11) 

1 .PP0F-02 

4.400E-03 

4.901E-03 

116.2 

12 

P(10) 

PI 12 ) 

1 . 88GF- 02 

4. 400C-03 

4.901E-03 

296. 2 

13 

Pill ) 

P( 13) 

1 .880E-02 

4.400E-03 

7.1 01E-03 

116.3 

14 

P( 12) 

P( 14) 

1.88CE-02 ' 

4.4C0F-03 

7. 101E-03 

296.3 

15 

P( 13) 

P ( 15 ) 

1.880E-02 

4.400E-03 

8.0C7E-03 

116. 3 

16 

P( 14) 

P( 16) 

1 .800F-02 

4.400E-03 

8.067E-03 

296. 3 

17 

P ( 6 1 ) 

P ( 1 7 ) 

2.950E-02 

3.800F-03 

6.238E-03 

288.6 

18 

P ( 6 1 3 

P(16) 

2.990E-02 

3. 800E-03 

6.2 33E-0 3 

108.6 

19 

P { 17 ) 

“ P( 19) 

2. 99 OF- 02 

3.800E-03 

5.32CE-03 

288.6 

20 

P 1 1 8 ) 

P ( 20 ) 

2.9SCF-02 

3. 800E-03 

5. 32BE-03 

100.6 

21 

P( 15) 

P (21 ) 

l .79CE-02 

4. 1 00F-03 

7.722L-03 

116. 3 

22 

P( 19) 

P ( 22 ) 

3.270E-02 

4. 1 00F-03 

3.640E-0 3 

208.6 

23 

P ( 20 1 

P ( 23 ) 

3.2 70E-02 

4. 100F-C3 

3.640E-03 

100.6 

24 

P 1 16) 

P ( 24 ) 

l .790F-02 

4. 100E-03 

7.722E-03 

296.3 

25 

P ( 2 l ) 

P ( 25 ) 

6 . 020E-01 

1. 700E-03 

3. 296E-03 

115.4 

26 

P ( 22 ) 

P( 29 ) 

6.020E-01 

l. 700E-03 

1 .898E-03 

295.0 

27 

P ( 23 ) 

P( 33 ) 

6.020E-01 

1.700E-03 

1. 898E-03 

1)5.0 

28 

P( 24) 

P ( 37 ) 

6. 020E-01 

1. 700F-03 

3.296C-03 

295.4 

29 

P ( 2 1 ) 

P ( 4 l ) 

6.360E-01 

0.0 

1.020E-03 

116.2 

30 

P ( 22 ) 

P (44) 

6.360E— 01 

0.0 

1.817E-03 

116.0 

31 

P C 23 ) 

P (47 ) 

6.360F-01 

0. 0 

1 .817E-03 

298.0 

32 

P ( 2 4 ) 

P ( 50 ) 

6.360E-01 

0.0 

1.020E-03 

296.2 

33 

P ( 2 l ) 

P( 53 ) 

2 .500F+00 

4. 500E-03 

5. 862E-04 

121.9 

34 

P 122) 

P( 55) 

2.500E+00 

4. 500E-03 

5.141E-04 

98.6 

35 

P( 23) 

P( 57) 

2.5CCF+00 

4. 500E-03 

5.141E-04 

278.6 

36 

P ( 24 ) 

P ( 59 ) 

2.500E +00 

4. 500E-03 

5.861E-04 

301.9 

37 

P ( 40 ) 

P ( 25 ) 

0. 0 

1. 550E— 02 

l. 143E-03 

104. 1 

38 

P ( 2 5 ) 

P ( 2 6 ) 

0.0 

1. 550E-02 

1 .183E-03 

106.9 

39 

P ( 26 ) 

P ( 27 ) 

0.0 

1. 550E-02 

1.635E-03 

208.4 

40 

P ( 27 ) 

P ( 28 ) 

0.0 

l . 550F-02 

4.905E-03 

294.1 

41 

P ( 28 ) 

P ( 2 9 ) 

0.0 

1. 550E-02 

1.660E-03 

284.5 

42 

P ( 29 ) 

POO) 

0.0 

1. 550E- 02 

6.663E-04 

288.2 

43 

POO) 

P (3 1 ) 

0.0 

l. 550E-02 

1.915E-03 

105.3 

44 

P ( 3 1 ) 

P( 32) 

0.0 

1.550E-02 

5.039E-03 

112.7 

45 

P ( 32 ) 

P ( 33 ) 

0.0 

l. 550F-02 

1.757E-03 

101. 0 

46 

P ( 33 ) 

P ( 34 ) 

0.0 

1.550E-02 

7. 608E-04 

99.6 

47 

P ( 34 ) 

P( 35) 

0.0 

1. 550E-02 

1.768E-03 

20B.8 


Figure H-3, Input and Output Data for Flowrates 
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48 

P ( 35 ) 

P ( 36 ) 

0.0 

1. 550E-02 

4.771E-03 

294.2 

49 

P ( 36 1 

P ( 3 7 ) 

0.0 

1. 550F-02 

1.2416-03 

279.2 

50 

P ( 3 7 ) 

P ( 38 ) 

0.0 

1. 550F-02 

1.2 77E-03 

281.9 

51 

P ( 38 ) 

P ( 39 ) 

0.0 

1. 550E-02 

1.4 956-03 

112.9 

52 

P ( 39 ) 

P ( 40 ) 

0.0 

1. 550F-02 

4. 6416-03 

115.8 

53 

P ( 52 ) 

P (41 ) 

0.0 

1. 140F-02 

7.91 5E-03 

113.2 

54 

P ( 4 1 ) 

P ( 42 ) 

0. 0 

1. 140E-02 

2 .5696-0 3 

106. 3 

55 

P< 42 ) 

P ( 43 ) 

0.0 

1. 140E-02 

2.506E-03 

283.9 

56 

P ( 4 3 ) 

P (44) 

0.0 

l. 140E-02 

8.0386-03 

292. 3 

57 

P ( 44 ) 

P( 45 ) 

0.0 

1. 140F-02 

' 1.659E-03 

287.1 

58 

P ( 4 5 ) 

P (46) 

0.0 

1. 140E-02 

' 2.9606-03 

100.9 

59 

P { 4 6 ) 

P (47 ) 

0.0 

1. 1406-02 

9. 0586-03 

110.9 

60 

P ( 47 ) 

P (48 ) 

0.0 

1. 14 OF — 02 

1.9G6E-03 

101 .0 

61 

P ( 4 8 ) 

P (49 ) 

0.0 

1. 14 0E-02 

2.5066-03 

283.9 

62 

P 149) 

P ( 50 ) 

0.0 

1. 140E-02 

B.132E-03 

291.6 

63 

P ( 5 3 ) 

P ( 5 1 ) 

0.0 

1. I 4 OF- 02 

2.8366-03 

28 2.3 

64 

P < 5 1 ) 

P ( 5 2 ) 

0.0 

1. 140E-02 

2.062E-03 

108.2 

65 

P 160) 

P< 53) 

0.0 

? . 900E-02 

3. 222E-03 

114.9 

66 

P( 53) 

P ( 54 ) 

0.0 

2. 900F-02 

9.475E-05 

143. 7 

67 

P ( 54 ) 

P ( 55) 

0.0 

2.900E-02 

3.5 32E-0 3 

293.0 

66 

P ( 55 ) 

P ( 5 6 ) 

0.0 

2.9 COE— 02 

2 . 7 5 9 E - 0 4 

79.6 

69 

P< 56 1 

P (57) 

0.0 

2. 900E-02 

3.668E-03 

111.2 

70 

P( 57) 

P ( 50 ) 

0.0 

2. 9006-02 

1.034E-04 

272.2 

71 

P ( 5 8 ) 

P ( 59) 

0.0 

2.9C0E-02 

3.354E-03 

292.8 

72 

P ( 59 ) 

P ( 6 0 ) 

0.0 

2. 900E-02 

2. 2 53E-04 

275.0 

73 

P ( 2 5 ) 

P IN ( l) 

l .54CE + 02 

2. 3706-02 

2.740E-03 

114. 4 

74 

P ( 26 ) 

P I N ( 2) 

1.540E+02 

2.3706-02 

2.242E-03 

106.4 

75 

P( 27 ) 

P IN ( 3) 

l .5408*02 

2.370E-02 

2.929E-03 

117. 1 

76 

P ( 20 ) 

P IN ( 4) 

l . 540E *02 

2.3 70E-02 

2 .957E-03 

299.8 

77 

P ( 29 ) 

P I N ( 5) 

1.54GE+02 

2. 3706-02 

2.202E-03 

291.1 

78 

P( 30 ) 

P IN ( 6) 

1. 540E*02 

2. 370E-02 

1 .935E-03 

285.2 

79 

P ( 3 1 ) 

P IN ( 7) 

1 . 540E *02 

2.370E-02 

2.769E-03 

298.0 

80 

P< 32) 

P IN ( 0) 

1. 540F*02 

2.37CE-C 2 

3.023E-03 

119.2 

81 

P< 33) 

P IN ( 9) 

1.540F+02 

2. 370E-02 

2.202F-03 

111.1 

82 

P( 34) 

P IN ( 10) 

l. 54CE+C? 

2. 370E-02 

1 .863E-03 

105.6 

83 

P ( 3 5 ) 

P IN ( 11) 

1.540E+02 

2. 370E-02 

2.609E-03 

1 19.0 

84 

P (36) 

P I N ( 12) 

I. 540E*02 

2 • 370E-02 

3. 332E-03 

298.7 

05 

P { 37 > 

PIN( 13) 

1.5406*02 

2. 37CE-02 

2.740E-03 

294.4 

86 

P ( 38 ) 

P IN ( 14) 

1 • 54CF *C2 

2 . 370E-02 

2. 171 F- 03 

286.8 

87 

P ( 39 ) 

P IN ( .15) 

1.5406*02. 

2.370E-02 

2.769E-03 

298.0 

88 

P ( 40 ) 

P IN i 16) 

1 . 54CF+02 

2.370F-02 

3.266F-03 

1 19. 2 

89 

P ( 4 1 ) 

P IN ( 17) 

7. 570E+01 

l. 170E-02 

5.860E-03 

116.5 

90 

P< 42) 

P IN ( 18) 

7. 57CE *01 

I . 17CE-02 

4.222E-03 

103.7 

91 

P ( 4 3 ) 

P IN ( 19 ) 

7. 5 70E + 01 

l. 170E-02 

5 .842E-03 

115. 5 

92 

P( 44) 

P IN (20) 

7. 5 7CE *Cl 

1. 17CE-02 

4. 6426-03 

292.6 

93 

P ( 4 5 ) 

DIN (21) 

7.5706*01 

1. 170F-02 

3.&44F-03 

281.7 

94 

P ( 4 6 ) 

P IN ( 22) 

7.5706*01 

l. 170E-02 

5.608E-03 

295.9 

95 

P ( 47 ) 

P INI 23) 

7.5706*01 

I . 17CE-02 

4.642E-03 

112.6 


Figure H-3. (Continued) 
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96 

P (48) 

P IN f 24) 

7.570E+01 

1. 170E-02 

3.4 06E-03 

101.5 

97 

P ( 49 ) 

P IN (2 5) 

7 . 570E+CI 

1. 17 CF -02 

5.014E-03 

115.9 

9B 

P ( 50 1 

P IN ( 26) 

7.570E+01 

1 . 170F-02 

5.860E-03 

29 6. 5 

99 

P ( 5 1 ) 

P IN (27) 

7. 570C *01 

l. 170E-02 

3.9H4E-03 

283.6 

100 

PI 52) 

P IN ( 28) 

7.570E «-0l 

l. 170F-02 

5.248E-0 3 

295. 4 

101 

P ( 53 ) 

P INI 29) 

1.280E+02 

l . 9SCE-02 

3. 37 5 E -03 

115.2 

102 

P ( 54 ) 

PIN (30) 

1.280E+02 

1.990F-02 

3.334E-03 

113.4 

. 103 

P ( 55 ) 

PIN! 31 ) 

1.260EVC2 

1 . 990E-02 

2.8 18E-03 

293.2 

104 

P ( 56 ) 

PIN I 32) 

1.280F+02 

1.990E-02 

3.1 32E-03 

29 3.6 

105 

P ( 57 ) 

P INI 33) 

1.2 80E* 02 

1.990E-02 

2.818E-03 

113.2 

106 

P I 58 ) 

P INI 34) 

1.2 8CE + 02 

l . 990E-02 

2.93 IE-0 3 

113. 9 

107 

P ( 59 ) 

PIN! 35) 

1.280EO2 

l . 990F-02 

3.375E-03 

29 5. 2 

108 

P ( 60 ) 

P INI 36) 

1 .2 3CE + C2 

1.990E-02 

3. 132E-03 

293.6 

XRL 

INJECTOR 






OXIDIZER SIDE, NO 

RING DAMS 





ANOVQIOUS SECOND T ANGENT I AL MODE 





• 




LB/SEC/PSI 

t FLOW/?' PI 

TOTAL 

INJECTOR FLOW 

~ 



1.2464E-01 

9. 1079 F— 01 

TOTAL 

VECTOR INJECTOR FLOW= 



1.24 16E-C1 

9.072SE-01 

TOTAL 

INJECTOR FLOW 

PROPORT IONE.D 

BY PC AMPLITUDES= 

1.2464E-C1 

9.1079 E— 0 1 

TOTAL 

VECTOR INJECTOR FLOW PROPORTIONED BY PC 

AMPLITUDES= 

1.24 16 E— 01 

9. 0 729 E— 01 


Figure H-3. (Continued) 


H-7 


PRES SUP. 

c VOLUME 

AC. VEl . 

MAGN ITUnE 

PHASE 

FLOWS IN 

FLOWS l 

OUT 

M JOE 

CU IN 

IN/SEC 

P S I /P S 1 

DEGREFS 






1 

4.040E+00 

5.3106+04 

2. 853E-01 

206.3 

1 


3 



2 

4. 04CE+00 

5.3L0E+04 

2.853E-01 

26. 3 

? 


4 



3 

4.04CE+00 

5. JlCE + 04 

5. 2436-01 

206. 3 

3 


5 



4 

4. 04 OF *-00 

5.3106+04 

5.243E-C1 

26.3 

4 


6 



5 

4.04CE+00 

5.31CE+04 

6.78CE-01 

2 06.3 

5 


7 



6 

4.040E+00 

5 .3 ICC+04 

6. 78CF-CI 

26.3 

6 


8 



7 

4. 140E+00 

5.3101+04 

7.2UE-01 

206.3 

7 


9 



8 

4.04CF+00 

5.31CE+04 

7.216F-CI 

26.3 

8 


10 



9 

4.04CE+00 

5. 3 ICE +04 

6.47SF-CI 

206. 3 

9 


1 1 



10 

4 .04 OF +00 

5.3106+04 

6. 4796-01 

26.3 

10 


12 



1 1 

4.040E+00 

5. 310E+04 

4. 688E-0A 

206.3 

11 


13 



12 

4.04CE+00 

5.310C+04 

4.688E-01 • 

26.3 

12 


14 



1 3 

4 < 940E + 0 0 

5.31 CE+C4 

2. 136E-01 

2 06.4 

13 


15 



14 

4.04 OE+OO 

5. 3106+04 

2. 136E-01 

26.4 

14 


16 



15 

4. 04 OE+OO 

5. 31. OF +04 

7.631E-02 

25. 8 

15 


21 



16 

4.040E *-00 

5.3106+04 

7.632E-02 

205. 8 

16 


24 



17 

4.200E+00 

5 . 3 10F+04 

l. 93 fE-Cl 

198.5 

17 


19 



18 

4. 200E+00 

5.310E+04 

1.936E-01 

18. 5 

18 


20 



19 

4 . 200E+00 

5.3 106+04 

3. 590E-01 

198.5 

19 


22 



20 

4. 2 CCE <-00 

5 . 3 1CE + 04 

3.59CE-01 

18. 5 

20 


23 



21 

7. 53CE+00 

5.310E+04 

3. 349E-01 

26.2 

21 


25 

29 

33 

22 

7.530E+00 

5.3106+04 

4.809F-01 

198.5 

22 


26 

30 

34 

23 

7. 53CE+00 

5.3 10E+0 4 

4. COSE— 0 1 

18.5 

23 


27 

31 

35 

24 

7.53GE+00 

5.3 10E+04 

3. 34SE-01 

206.2 

24 


28 

32 

36 

25 

1.190E+00 

5 . 3 10E+04 

3. 807E-01 

25. 8 

25 

37 

38 

73 


26 

9. 770 E— 01 

5.310E+04 

5. 292F-01 

23.3 

38 


39 

74 


27 

9.770E-0 1 

5. 3 ICE + 04 

3. 235E-01 

26.4 

39 


40 

75 


28 

9.770E-01 

5.3 1 CE + 04 

2.981E-01 

201. 8 

40 


41 

76 


29 

l.l 90E+OC 

5.3 1 OE+04 

5.C72E-01 

198. 7 

41 

26 

42 

77 


30 

9. 770 E— 01 

5.3 10E+04 

5.916E-01 

198. 7 

42 


43 

78 


31 

9.770E-01 

5.3 10E + 04 

3.49FE-01 

201. 0 

43 


44 

79 


32 

9.77CE-01 

5 . 31 CE+04 

2. 8U5E-01 

24. 7 

44 


4 5 

80 


33 

1 • 190 E+00 

5 . 3 10E+04 

5. 072F-01 

18.7 

45 

27 

46 

81 


34 

9.770E-01 

5. 3 10E+04 

6.025F-01 

17.3 

46 


47 

82 


35 

9.770E-0 l 

5. 310F+04 

3. 788E-C1 

16.4 

47 


48 

83 


36 

9.770E-0 l 

5.3 10E+04 

2. 344E-01 

216. 8 

48 


49 

84 


37 

1. 190E+00 

5. 31 06 +04 

3. 807E-01 

205.8 

49 

28 

50 

85 


30 

9. 770 E- 01 

5.3106+04 

5.390F-01 

201.6 

50 


51 

86 


39 

9.7706-01 

5.3 10E+04 

3.4S8E-01 

201. 0 

51 


52 

87 


40 

9.770E-01 

5. 310E+04 

2. 407E-01 

32. 7 

52 


37 

88 


41 

1.37 CE+OO 

5.3 10E +04 

3.34SE-01 

26. 0 

53 

29 

54 

89 



Figure H-4, Input and Output Data for Pressure Nodes 
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42 

1.37 CE + OO 

5.3 10E+04 

5. 74CE-01 

22. 0 

54 


55 

90 

43 

1.370E+00 

5.31 OE *04 

3.446E-01 

2 7.5 

55 


56 

91 

44 

1 .370E + 00 

5.310E*04 

4. 6C7E-01 

193. 7 

56 

30 

57 

92 

45 

1.37CE+0C 

5.3 10F+04 

6. 351E-01 

198. 3 

57 


58 

93 

46 

1 . 370E+00 

5.310F*04 

3. 636E-01 

203.9 

58 


59 

94 

47 

1 . 370E + 00 

5. 310F+04 

4.807E-01 

18. 7 

59 

31 

60 

95 

48 

1.370E+00 

5 . 3 ICE *04 

6.571F-01 

16. 6 

60 


61 

96 

49 

1 .3 7CE *00 

5.3 10E + 0 4 

4.241E-01 

18.1 

61 


62 

97 

50 

l . 3 70F +CC 

5. 310 F. *04 

3.34 9F-01 

206.0 

62 

32 

63 

90 

51 

1.37CE+00 

5.310E+04 

5. 94QE-01 

200.0 

63 


64 

99 

52 

1.370E+00 

5. 3 I DF *04 

4.029E-01 

200. 8 

64 


53 

100 

53 

8.660E-0 1 

5.3106+04 

3.565E-01 

26.3 

65 

33 

66 

101 

54 

6.660E-01 

5.310F+04 

3. 766E-01 

27. 8 

66 


67 

102 

55 

8.66CE-01 

5.310E+04 

4.626E-01 

199. 1 

67 

34 

68 

103 

56 

6.660E-0 1 

5.310F+04 

4.069E— 01 

203.6 

68 


69 

104 

57 

8.66CE-01 

5.310E+04 

4. 62CE-01 

19. 1 

69 ’ 

35 

70 

105 

58 

6. 660 E- 01 
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Figure H-4 . (Continued) 


17 IP 


K-9/H-10 



